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Abstract 

We discuss the minimal supersymmetric U(1)b— l x U(1)r extension of the standard model. 
Gauge couplings unify as in the MSSM, even if the scale of U(1)b-l x U(1)r breaking is as 
low as order TeV and the model can be embedded into an SO(W) grand unified theory. The 
phenomenology of the model differs in some important aspects from the MSSM, leading potentially 
to rich phenomenology at the LHC. It predicts more light Higgs states and the mostly left CP-even 
Higgs has a mass reaching easily 125 GeV, with no constraints on the SUSY spectrum. Right 
sneutrinos can be the lightest supersymmetric particle, changing all dark matter constraints on 
SUSY parameter space. The model has seven neutralinos and squark/gluino decay chains involve 
more complicated cascades than in the MSSM. We also discuss briefly low-energy and accelerator 
constraints on the model, where the most important limits come from recent Z 1 searches at the 
LHC and upper limits on lepton flavour violation. 
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I. INTRODUCTION 



Within the minimal supersymmetric extension of the standard model (MSSM) the gauge 
couplings unify nearly perfectly around an energy scale of approximately ~ 2 x 10 16 GeV, 
if SUSY particles exist with masses of the order of 0(1) TeV. Extending the MSSM with 
non-singlet superfields tends to destroy this attractive feature, unless (a) the additional fields 
come in complete SU (5) multiplets or (b) the standard model gauge group is extended too. 
Here we study a model in which the SM group is enlarged to SU (3) c x SU(2) l x U(1)b-l x 



U(V)r. It is a variant of the models first proposed in |1[ and later discussed in more detail 
in @. 

Our main motivation for studying this model can be summarized as: (i) It unifies, in the 
same way the MSSM does, even if the scale of U(1)b-l x U(1)r breaking is as low as the 
electro-weak scale; (ii) it can be easily embedded into an SO (10) grand unified theory; (iii) 
it has the right ingredients to explain neutrino masses (and angles) by either an inverse [3J] 
or a linear [4,0] seesaw; (iv) it allows for Higgs masses significantly larger than the MSSM 
without the need for a very heavy SUSY spectrum |s| and (v) it potentially leads to rich 
phenomenology at the LHC. 

With the data accumulated in 2011 both ATLAS 0] and CMS [s| have seen some in- 
dications for a Higgs boson with a mass of roughly ~ 125 GeV. This result, perhaps 
unsurprisingly, has triggered an avalanche of papers studying the impact of such a relatively 
hefty Higgs on the supersymmetric parameter space 0-0. The general consensus seems to 
be, that the MSSM can generate ~ 125 GeV only if squarks and gluinos have masses in 
the multi-TeV range. While this is, of course, perfectly consistent with the lower bounds on 
SUSY masses obtained from Jfix searches at the LHC [4J, |45(, such a heavy spectrum could 
make it quite difficult indeed for the LHC to find direct signals for SUSY. 

There are, of course, several possibilities to circumvent this conclusion. First of all, it is 
well-known that the loop corrections to h° are dominated by the top quark-squark loops. 
Thus, little or no constraints on sleptons and on squarks of the first two generations can in 
fact be derived from Higgs mass measurements, once the assumption of universal boundary 
conditions for the soft SUSY parameters is abandoned. Second, in the next-to-minimal SSM 
(NMSSM) the h° can be heavier than in the MSSM due to the presence of new F-terms 
from the additional singlet Higgs 18|, |43j , especially in models with non-universal boundary 
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27 . And, 



conditions for the (soft) Higgs mass terms 25J or in the generalized NMSSM 
third, in models with an extended gauge group additional D-terms contribute to the Higgs 
mass matrices, relaxing the MSSM upper limit considerably 46m0| . This latter possibility is 
the case we have studied in a previous paper jsj using the minimal U(1)b-l x U(1)r models 
of [l| . Here, we extend the analysis of p] , including both Higgs and SUSY phenomenology. 

Due to the extended gauge structure the model necessarily has more Higgses than the 
MSSM. Near D-flatness of the U(1)b-l x U(1)r breaking then results in one additional 



light Higgs, h BLR [6J]. Mixing between the MSSM h° = h° L and h BLR enhances the mass 
of the mostly MSSM Higgs and, potentially, affects its decays. This is reminiscent to the 
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situation in the NMSSM, where an additional light and mostly singlet Higgs state seems to 



be preferred [18|, |43| if the signals found by ATLAS [7| and CMS [8| are indeed due to a 125 
GeV Higgs. 

The MSSM-like h° L in our model can have some exotic decays. For example, the h° L will 
decay to two lighter Higgses, if kinematically possible, although this decay can never be 
dominant due to constraints coming from LEP. The model also includes right-handed neu- 
trinos with electro-weak scale masses and there is a small but interesting part in parameter 
space where m z o < m VR < m fe o, where the Higgs decays to two neutrinos. These decays 
always lead to one light and one heavy neutrino, with the latter decaying promptly to either 
W ± l T or Z°u. (Mostly right) sneutrinos can be lighter than the h° L , in which case the Higgs 
can have invisible decays. 

The SUSY spectrum of the model is also richer than the MSSM: It has seven neutralinos 
and nine sneutrino states. These additional sneutrinos can easily be the lightest supersym- 
metric particle (LSP) and thus change all the constraints on SUSY parameter space, usually 
derived from the requirement that the neutralino be a good dark matter candidate with the 



correct relic density 5JJ 1 . Even though the lightest sneutrino can also be the LSP in the 



MSSM, direct detection experiments have ruled out this possibility a long time ago [52 
In SUSY decays, within the MSSM right squarks decay directly to the bino-like neutralino, 
leading to the standard missing momentum signature of supersymmetry. Due to the ex- 
tended gauge group, right squarks can decay also to heavier neutralinos, leading to longer 
decay chains and potentially to multiple lepton edges 2 . Decays of the heavier neutralinos 
also produce Higgses, both the h° L and the h° BLR appear, with ratios depending on the right 
higgsino content of the neutralinos in the decay chains. 

The rest of this paper is organized as follows. In the next section we discuss the setup of 
the model, its particle content, superpotential and soft terms and the symmetry breaking. 
The phenomenologically most interesting mass matrices of the spectrum are given in section 
II III where we also discuss numerical results on the SUSY and Higgs mass eigenstates. Here, 
we focus on Higgs and slepton/sneutrino masses, which are the phenomenologically most 
interesting. In section |IV] we define some benchmark points for the model and discuss their 
phenomenologically most interesting decay chains. We then close with a short summary. 
In the appendix we give mass matrices not presented in the text, formulas for the 1-loop 
corrections in the Higgs sector and more information about the calculation of the RGEs, 
including anomalous dimensions as well as the 1-loop functions for gauge couplings and 
gauginos. 



1 Also for the case of a neutralino LSP, constraints on the SUSY parameter space from dark matter can 
change in case the right Higgsino is light. 



Longer SUSY cascades from larger number of neutralino states have also been discussed in 53, [54 [ . 
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Superfield 


SU(3) C x SU(2) L x U(1) R x U{1) B -l 


Generations 




Q 

d c 

u c 
L 
e c 
u c 
S 


(3,2,0,+i) 
(3, 1, +2, — g) 

( 3 . 1 . 7) « ~F 1 

V * * 2 ' 6 / 

(1,2,0,-|) 
(1,1,0,0) 


3 
3 
3 
3 
3 
3 
3 




H u 
H d 

XR 
XR 


(1,2, +|,0) 
(1,2,-1,0) 

(i,i,+i-i) 
(M, -§.+§) 


1 
1 
1 
1 



TABLE I: The Matter and Higgs sector field content of the U(1)r x U(1)b-l model. Generation 
indices have been suppressed. The S superfields are included to generate neutrino masses via the 
inverse seesaw mechanism. Under matter parity, the matter fields are odd while the Higgses are 
even. 

II. THE MODEL: SU{3) C x SU(2) L x U(1) B -L x U(1) R 

In this section we present the particle content of the model, its superpotential and discuss 
the symmetry breaking. We consider the simplest model based on the gauge group SU(3) C x 
SU(2) L x U(1) R x U(1) B -L- We will call this the mBLR model below. As has been shown 
in [l| it can emerge as the low-energy limit of a certain class of SO(10) GUTs broken along 
the "minimal" left-right symmetric chain 

5O(10) SU(3) C x SU(2) L x SU(2) R x U(1) B -l (1) 
-> SU(3) C x SU(2) L x U(1) R x U{1) B -l. 

The main virtue of this setting is that an MSSM-like gauge coupling unification is achieved 
with a sliding U(l) R x U(1)b-l breaking scale, i.e. this last stage can stretch down even 
to the electro- weak scale. Different from the previous works [H, 0] , we assume that the first 
two breaking steps down to U(1) R x U(1)b-l happen both at (or sufficiently close to) the 
GUT scale. This assumption is used only for simplifying our setup, it does not lead to any 
interesting changes in phenomenology. 

A. Particle content, superpotential and soft terms 

The transformation properties of all matter and Higgs superfields of the model are sum- 
marized in table HI Apart from the MSSM fields, in the matter sector we have v c and S. The 
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former are necessary in the extended gauge group for anomaly cancellation, 3 while the fields 
S are included to explain neutrino masses by either an inverse [2I or a linear 0, \^ seesaw 
mechanism. Our Higgs sector, including the new fields xr and xr-, is the minimal one for 
the breaking of U(1)b-l x U{\)r to U{1)em- 

The fields xr an d Xr can be viewed as the (electric charge neutral) remnants of SU(2)r 
doublets, which remain light in the spectrum when the SU(2)r gauge factor is broken by 
the vev of a B — L neutral triplet down to the U(1)r The presence of xr an d Xr makes 
it necessary to introduce an extra Z^ 1 matter parity, since otherwise .R-parity is broken in 
a potentially disastrous way, once these scalars acquire vacuum expectation values. This 
Z^ 1 is not a particular feature of our setup; it is always needed in models where U(1)b-l is 



broken with doublets 55 



The relevant -R-parity and Z^ 1 conserving superpotential is given by 

W = Wmssm + W s . (2) 

Here, 

Wmssm = Y u u-QH U - Y d d c QH d - Y e e c LH d + /iH u H d (3) 
W s = Y u u c LH u + Y s u c xrS - hrXrXr + l^sSS. 

where Y e , Y d and Y u are the usual MSSM Yukawa couplings for the charged leptons and the 
quarks. In addition there are the neutrino Yukawa couplings Y v and Y s ; the latter mixes the 
v c fields with the S fields giving rise to heavy SM-singlet pseudo-Dirac mass eigenstates. The 
term fi R is completely analogous to the MSSM fi term. Note that the term fi S is included 
to generate non-zero neutrino mass with an inverse seesaw mechanism. However, as always 
is done in inverse seesaw, we assume that fis is much smaller than all other dimensionful 
parameters of the model. Apart from neutrino masses themselves it will therefore not affect 
any of the mass matrices (or decays) of our interest. 

Note that, besides the role it plays in neutrino physics, the Y s coupling is relevant also 
for the Higgs phenomenology at the loop level as it enters the mixing of xr an d Xr Higgs 
fields with the SU(2)l Higgs doublets as well as the RGEs for xr, see below. 



Following the notation and conventions of 56[ the soft SUSY breaking Lagrangian reads 



Vsoft = Yl Ktth + ( M - X « X - + T ^rQ H u - T d d* R QH d + T v v* R LE u 

ij a 

-T e e* R LH d + B»H u H d - B^xrXr + T s D rX rS + B^SS + h.c^j . (4) 

The first sum contains the scalar masses squared and the second sum runs over all gauginos 
for the different gauge groups (called \rl, Ar, \ l L and Xq in the following) and the second one 
contains the scalar masses squared. While B^ s is in principle a free parameter, a naive order 



v c is automatically part xithe theory due to its 50(10) origin. 
In the normalization of 



55| doublets have U(1)b-l = b i-e are "odd" under B-L. 
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of magnitude expectation for it is B^ s ~ Hs m susY- Thus, one expects that B^ s is much 
smaller than all other soft terms and can be safely neglected, see discussion of sneutrinos 
below. 

To reduce the number of free parameters, in our numerical studies we will consider a 
scenario motivated by minimal supergravity. This means that we assume a GUT unification 
of all soft-breaking sfermion masses as well as a unification of all gaugino mass parameters 

mlSij = m 2 D Sij = m^Sij = m 2 Q 5ij = m 2 ^ = m|5 i3 - = m^^- (5) 

Mi /2 = M BL = M R = M 2 = M 3 

Also, for the trilinear soft-breaking coupling, the ordinary mSugra conditions are assumed 

Ti = A Yi, i = e,d,u,v,s . (6) 

The GUT scale is chosen as the unification scale of gsL, 9r and gL, while we allow g 3 to be 
slightly different, exactly as in the MSSM. A complete unification is assumed to happen due 
to GUT threshold corrections. For the remaining soft parameters in the Higgs sector, m 2 H , 
m 2 Hu1 m| and and B m , we have implemented two different options. These 

are discussed in section III Bt 

The presence of two Abelian groups gives rise to gauge kinetic mixing 



This is allowed by gauge and Lorentz invariance 57J], as F B i,AU/ and F R, ^ U are gauge 



invariant, see e.g. [58]. Even if U(1)r and U(1)b-l are orthogonal in SO (10) the kinetic 
mixing term will be induced during the RGE running below the SU(2)r breaking scale 
because the light fields remaining below the GUT scale can't be arranged in complete 50(10) 
multiplets: while all matter fields form three generations of 16-plets, xr and \r induce off- 
diagonal elements already in the 1-loop matrix of the anomalous dimensions defined by 
Irbl = j^TrQ R Q B _ L . The matrix reads 

3 




N = diag(l, w|) contains the GUT normalization of the two Abelian gauge groups. Our 



implementation follows the description of 59], where it is shown that terms of the form as 



in eq. (|7|) can be absorbed in the covariant derivative by a re-definition of the gauge fields. 
Therefore, we are going to work in the following with covariant derivatives of the form 

D li = d li -iQ%GA li , (9) 

where Q$ is a vector containing the charges of the field $ with respect to the two Abelian 
gauge groups and G is the gauge coupling matrix 

G= f 9R 9RBL \ (1Q) 
\ QBLR gBL I 
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contains the gauge bosons = (A R , A^ L ) T . Since the off-diagonal elements in eq. (jHJ) 
are negative and roughly one order smaller than the diagonal ones, it can be expected that 
the off-diagonal gauge couplings at the SUSY scale are positive but also much smaller than 
the diagonal ones. This is in some contrast to models in which kinetic mixing arises due to 
the presence of U(l)y X U(1)b-l 60] • in addition, a mixing term of the form 

M BLR X BL X R (11) 

between the two gaugino Xrl and Xr will be present [6l|. Since we have chosen the SU(2)r 
breaking scale to be very close to the GUT scale we demand as additional boundary condi- 
tions that the new parameters arising from kinetic mixing vanish at the GUT scale, i. e. 

9RBL = 9BLR = , M B LR = . (12) 

For more details on UJV) mixing and its physical impact we refer the interested reader also 



60 



62-64 . Our focus will be on the additional terms in the scalar mass 



to recent papers 
matrices due to the presence of non-diagonal couplings. 



B. Tadpole equations and boundary conditions 

The U(1)r x U (1) b-l gauge symmetry is spontaneously broken to the hypercharge U(1)y 
by the vevs v XR and v XR of the scalar components of the \R an d Xr superfields while the 
SU "(2) l®V (l)y —> U(1)q is governed by the vevs v d and v u of the neutral scalar components 
of the SU{2) L Higgs doublets H d and H u up to gauge kinetic mixing effects. One can write 

XR = ^ ( Gr + i( ^R + v xr) > Xr = ~j= ipR + i<PR + v XR ) » ( 13 ) 

H° d = ^( a d + i<Pd + v d ) , H u = ^ iPu + ifu + v u ) , (14) 

where the generic symbols a and ip denote the CP-even and CP-odd components of the 
relevant fields, respectively. 

The minimum conditions for the four different vevs can be written at tree-level as 

\a lr M - vl) + l -A LR ^vl R - v 2 XR )j (15) 
\a L rM - vl) - \a lr ^vI r - (16) 

hn = ~ B m v XR + v xr ( m L + M 2 + 1 A lrA v Ir ~ v Ir) + \ A lrM - ^u)J ( 17 ) 

t XR = - B m v XR + v xr ( m L + M 2 - I a lr,i(v Xr - v 2 XR ) - ^A LRj2 (v 2 d - vl)J (18) 
where we defined 

Alr,i = 9bl + 9r + 9blr + 9rbl ~ ^9r9blr — ^9bl9rbl 

A 2,2 

A-lr,2 — 9r + 9rbl ~ 9r9blr — 9bl9rbl 

A L R,3 = 9l+9 2 R + 9 2 RBL- (19) 



td = 


-B^Vu + v d yn 2 Hd 4 


- H 




-B^v d + v u ( m 2 Hu H 


- \iA 
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For the vacuum expectation values we use the following parameterization: 



< = <n + 4n> J = vl + Vl (20) 

V V D ^ Vu 



XR XR 

j 

tan f3 R = — , tan (3 

V XR v d 



The tadpole equations can analytically be solved for either (i) (fi,B^,fi R , B^ R ) or (ii) 
(fi, Bp, m^, R , m^ R ) or (iii) (m 2 Hd , m 2 Hu , m^, fn^). Option (i) can be considered the minimal 
version. We call this option CmBLR (constrained mBLR), since it allows to define boundary 
conditions for all scalar soft masses, = m% n = rriQ and m^ fl = m£ = at mcuT, 
reducing the number of free parameters by four. This assumption, however, leads to some 
important constraints on the parameter space, as we will discuss next. Options (ii) and 
(iii) are more flexible. Option (ii) is similar to the CMSSM with non-universal soft masses 



(NUHM) 65M671]. albeit the non-universality is only in the B — L sector. We will call this 
the x.R m BLR (non- universal xr masses mBLR), and most of our numerical results are based 
on this option. We mention option (iii) for completeness, but we have not used it in our 
numerical studies. 

As will be shown in section IIII A\ the mass of the Z'-boson in the mBLR model is 
approximately given by 

1 



rrir 



-7 A LR,lV R (21) 



4 

We can use this expression and eqs (1I5]) - (ITB]) to obtain an approximate relation between mz> 



and [ir, m 2 , m 2 , and tan (3r. This leads to 



m % s -2( W » + <) + f ^cos(2«^|±i + A^^j (22) 
where Am 2 = m% — m 2 . We can roughly estimate Am 2 , if we make a mSugra-like 

XR XR XR XR 

assumption for the boundary conditions, m? H = fn\ R = wig at the GUT scale. The running 
value of Am^ R can then be found by a one-step integration of the RGEs at 1-loop level as: 



A>»L - j-pf-fttttf )(*>■? + AD log (23) 

with T s ~ A Y S . As eq. (!23j) shows, with these assumptions Am^ > and the condition 
that mz' of eq.( |22l) has to fulfill the experimental lower bound will define an excluded area 
in the 3-dimensional parameter space [Tr(Y s YJ), ta.n (3R,m 2 RGE ], where m RGE = (3mo + ^4q) . 
If we assume in addition that Y s is small enough to remain perturbative anywhere between 
the weak and the GUT scale, a lower bound on m 2 RGE as a function of tan (3r — 1 will result 
in the CmBLR. 

This can be understood in more details as follows. In the CmBLR Am^ s > 0, as shown 
by eq. (I2"3l) and the last term in eq. (|2"2"|) is positive only if tan (3r > 1. Since cos(2/3) < 

2 

the second term in eq. f[2~2"l) is positive only if tan (3r < 1. If Am^ B > \^f-v 2 cos(2/3)|, only 
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FIG. 1: Constraints on the CmBLR parameter space from the condition of correct symmetry 
breaking, to the left: v R = 5 TeV, to the right v R = 8 TeV. In both plots M 1/2 = 1000 GeV, 
tan/3 = 10 and Aq = 0. Just above the lines for fiR = —200 GeV \hr\ 2 = 0, i.e. larger values 
of tan/3ij do not lead to consistent solutions of the tadpole equations (for fixed mo and Aq). For 
detailed explanation see text. 



solutions with tan/?# > 1 can be found. Since finally |yUj?| 2 must be \[ir\ 2 > and m| > 
in the CmBLR we get the constraints on the parameter space shown in fig. [TJ Here we show 
for two choices of vr contour lines of fiR in the plane (tan [3 R ,m ). Just above the lines for 
fiR = —200 GeV |/!_r| 2 = 0, i.e. larger values of tan (3r do not lead to consistent solutions 
of the tadpole equations (for fixed m and A ). This restricts the model to values of tan/^ 
very close to 1, as is clearly demonstrated in the figure. Note that for low values of mo the 
constraints on the viable region of tan (3 R actually becomes stronger. 5 

No such constraint on ttlq and Aq exists in the x/? m BLR, since here Am^ is a free 
parameter. However, if (Am^/m 2 ^) <C 1, values of tan f3 R very close to 1 are preferred by 
eq. ([22]) in both, the CmBLR and the Xi?mBLR. 



III. MASSES 



In this section we give the most important mass matrices of the model at tree-level. In 
the numerical calculations we take also the 1-loop corrections jsij into account, see appendix 
for more details. The numerical implementation of the model has been done using SPheno 



69l. l70j , fo r which the necessary subroutines and input files were generated using the package 



SARAH 



71-73 



The used model files are included in the public version 3.1.0 of SARAH. 



5 tan/3i? ~ 1 is also needed for a spectrum without tachyons since the additional D-terms can give large 
negative contributions to the sfermion masses, see below. 
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A. Gauge bosons 



In the basis (W°, Bb-l, B r ) the mass matrix for the neutral gauge bosons reads at tree- 
level 

\ 



Mi 



vv 



1 

4 



g 2 y 



-9l9rblv 



9l9rv 



-9l9rblV 2 9rbl v2 + 9bl v r 9r9rblV 2 - 9r9blV r 



-9l9rv 2 9r9rblV 2 



9r9blv r 



9%v 2 + g 2 R v 2 R 



(24) 



/ 



where 
From eq. 



9bl = (gBL - 9rbl) , 9r = (9r ~ 9blr) ■ (25) 
the masses of the photon, the Z and the Z' can be calculated analytically 



in 



z,z> 



- I Av 2 + Bv 2 R T v 2 R < 



R 



R 



B 



(26) 



with 
B 



2,2,2 

9l + 9r + 9rbl 

9bl + 9r + 9 blr + 9rbl — 2>9blr9r — ^9rbl9bl 



C = g\{9R ~ 9blr) 2 + 9 2 bl(9l + 9r) ~ ^9bl{9\ + 9blr9r)9rbl + (g 



Expanding eq. [26] in powers of v 2 /v R , we find up to first order 



rrir 



Cv 2 
AB 



rrir 



(AB - C)v 2 + B 2 v 



2„,2 
R 



AB 



In the limit gBLR = and gRBL = we then get 



rrir 



(g 2 B L gl + 9bl9r + 9l9rY 



±(g 2 



BL 



g\) 



rrir 



gW 



BLR 



gl)g 



RBL 



(27) 



(28) 



<g 2 BL + g R ) 4 



\-Wbl + 9>1- (29) 



ATLAS has recently published updated lower limits on Z' searches 74j. Our Z' corresponds 
to the Z x in the notation of 75j, i.e. 74J gives a lower limit of Z' > 1.8 TeV, which 
corresponds to roughly v R > 5 TeV for our choice of couplings 6 , see, however, the discussion 
in section ITV CI 



B. Higgs bosons 

1. Pseudoscalar Higgs bosons 

At the tree level we find that in the (ifa, (p u , ft R , tp R ) basis the pseudoscalar sector has a 
block-diagonal form and reads in Landau gauge 



6 The condition that the gauge couplings reproduce correctly the standard model hypercharge, plus the 
assumption of unification lead to values of roughly gBL ~ 0.57, gn ~ 0.45, gBLR ~ 0.014 and gRBL ~ 0.012 
at the SUSY scale. 
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M AA 



M\ A r 



M 2 A A,R 



(30) 



with 



From these four states two are Goldstone bosons which become the longitudinal parts of 
the massive neutral vector bosons Z and a Z' . In the physical spectrum there are two 
pseudoscalars A and A° R with masses 

m 2 A = ^(tan/3 + l/tan/3) , m A R = B^ R (tan(3 R + I/tan/?/?) . (32) 



2. Scalar Higgs bosons 



The tree-level CP-even Higgs mass matrix in the (er^, a u , a R , ctr) basis reads 



m LR m RR 



where 



(si + I^IblK^ + -§ (m\ + (#f + ^| BL )w 2 ) s 2 p , 

' ii 1 A™2 , ( „2 , 1 „2 \„.2N „ ^„2 , 1 „2 \„.2„2 , ™2 „2 ' V / 




2 „2 
A c /3 



{qrQr- 9bl9rbl)vv r c I3 c I3r -(g R g R - gBLg R BL)VV R CpSp R . 

m LR — T /- ~ \ /- ~ \ ' l^^J 



-(9r9r - 9bl9rbl)w r s i3 c i3r (g R g R - gBLg R BL)vv R spSp R 



(36) 



s x = sin(x), c x = cos(x) (z = (3, (3 R ,2(3,2(3 R ), g% = {g\ + g 2 R )/A, ~g\ R = (g 2 BL + g 2 R )/A. 
The matrix m 2 LL contains the standard MSSM doublet mass matrix. To see this explicitly 
one has to integrate out the additional Higgs fields in the v R — > oo limit which yields a 
shift in the gauge couplings such that the MSSM limit is achieved. m RR corresponds to 
the U(1) R x U(1)b-l Higgs bosons and m 2 LR provides the essential mixing between the two 
sectors. 

Note that it is straightforward to show that the determinant of the mass matrix eq. ( 133]) 
goes to zero, whenever one of the parameters ((tan/3 — l),rriA, (tan/3^ — l),m^ H ) goes to 
zero. One can also calculate analytically that in the limit of v R — > oo the lightest eigenvalue 
of eq. fl3"3l obeys the MSSM tree- level limit for h°. For finite v R corrections to m 2 LL appear, 
of the order of g R v 3 /v R , which lead to a shift in the lightest eigenvalue. Thus the MSSM 
tree-level upper bound of m^f 6 < m z o for the lightest Higgs can be violated. 
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FIG. 2: Example plot for the masses (left) and "leftness" (right) of two lightest eigenvalues of the 
CP-even Higgs sector as a function of vr for fixed choices of the other parameters: ttiq = 250 GeV, 
Mi/2 = 800 GeV, tan = 10, A = 0, tan p R = 0.94, (j,r = -800 GeV, m AR = 2350 GeV. 

3. Numerical examples 

In fig. [2] we show the two lightest Higgs boson masses (to the left) together with 

Ki = Rl + Rl (37) 

(to the right) as a function of Vr. Here % = 1, 2 labels the light Higgs scalars in the model 
and .Rjj is the rotation matrix which diagonalizes the CP-even Higgs sector. Note that the 
quantity TZ 2 Li , which reaches one in the MSSM limit, is a rough measure of how much the 
corresponding Higgs with index % resembles an MSSM Higgs boson. We will call this leftness. 
Roughly speaking, the smaller this quantities is, the smaller is the z-th Higgs coupling to 
the Z- and VF-bosons, implying a reduced production cross sections at LEP, Tevatron and 
the LHC. 

Since the MSSM Higgs and the two additional Higgs xr an d Xr are both charged under 
U(1)r the two lightest Higgs states mix due to additional D-terms in the CP-even Higgs 
matrix (see eq. (|34"|) ). Thus, both the masses and the mixing of the two lightest Higgs hi 
and h 2 depend strongly on Vr, see fig. El In this example, up to approximately Vr—6 TeV 
h\ is mostly the singlet Higgs whereas h% is the MSSM-like Higgs. For larger vr a level- 
crossing occurs and the situation is reversed. Note that, although mo, M1/2 and Aq have 
rather moderate values in this example, hi has a mass of the order of rrih 2 — 125 GeV for 
vr ~ 5 — 6 TeV, i.e. the D-terms have shifted the MSSM-like Higgs mass into the region 
preferred by ATLAS and CMS. 

In fig. [3] masses and "leftness" of the two lightest Higgs eigenstates are plotted against 
tan/?R. For tan/^R close to tan/?R = 1 one gets a very light singlet Higgs as expected, see 
discussion above. As is the case when varying vr a level- crossing appears also when tan /3r 
is changed. In this figure the Xi? m BLR version was used, thus we can put tan/3/? < 1. Note, 
however, that the masses of hi and h 2 show a behavior which is symmetric with respect to 
I tan /3r — 1\. Again the figure demonstrates that the MSSM limit of the lightest Higgs mass 
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FIG. 3: Example plot for the masses (left) and "leftness" (right) of two lightest eigenvalues of 
the CP-even Higgs sector as a function of tan (3r for fixed choices of the other parameters: mo = 
250 GeV, M 1/2 = 800 GeV, tan/3 = 10, A = 0, v R = 6000 GeV, lir = -800 GeV, m AR = 
2350 GeV 
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FIG. 4: Example plot for the masses (left) and "leftness" (right) of two lightest eigenvalues of the 
CP-even Higgs sector as a function of lir for fixed choices of the other parameters: mo = 250 GeV, 
M1/2 = 800 GeV, tan/3 = 10, A = 0, VR = 6000 GeV, tan/3 R = 0.94, m Aa = 2350 GeV 



can be violated at the expense of a reduced coupling of the MSSM-like state to SM gauge 
bosons. 

1-loop corrections play in general an important role, not only for the MSSM-like Higgs but 
also in the singlet sector. This can be seen in fig. |4j where we plot masses and mixings of hi 
and hi versus the parameter llr. Increasing or decreasing lir, respectively, changes the mass 
of the mostly-singlet Higgs by considerable factors. In fact, for larger values of \lir\ one can 
get easily a negative mass squared for hi, which is of course forbidden phenomenologically. 
The importance of lir stems from 1-loop contributions to the Higgs mass matrix with a 
higgsino-right in the loop. Loop corrections for the mostly-singlet Higgs are, in fact, even 
more important numerically than for the MSSM-like Higgs and many points which are 
allowed at tree- level lead to tachyonic states, once 1-loop corrections are taken into account. 

Finally we note, that in the plots in this section we have not shown the regions excluded 
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by LEP or the LHC searches, since we were interested only in showing the parameter de- 
pendencies of our numerical results. In the study points of the next section, however, we 
have taken care that our points survive all known experimental constraints. 



C. Neutrinos 

The mBLR model contains beside the usual three left-handed neutrinos six additional 
states which are singlets with respect to the SM group. The corresponding mass matrix is 
in the basis (ul, z/ c , S) given by 



m. 



( 
\ 

This matrix is diagonalized by U u : 








— v Y 



±v Y T 









-^-v Y T 



MS 



(38) 



U v >*m„U 



m 



dia 



(39) 



Eigenvalues for the three light (and mostly left-handed) neutrinos can be found in the seesaw 
approximation as: 



off 



(40) 



Neutrino data implies that either Y v and/or \ig is small and in inverse seesaw the smallness 
of neutrino mass is attributed to the smallness of the latter. As we will discuss in section 
IIV Al the bounds on rare lepton decays imply that the off-diagonal terms of Y s and Y v have 
to be small compared to their diagonal entries, unless their diagonal values are small too. 

The smallness of us implies that the six heavy states form three "quasi-Dirac" pairs. For 
vanishing off-diagonal entries in Y s and Y v a good estimate of the masses of the heavy states 
is: 

(41) 



m u h ,ii 
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D. Sparticles 



1 . Neutralinos 

The mass matrix of the neutralinos reads in the basis (Xbl, X l , h% h^, Xr, Xr, Xr)- 



/ 


M BL 





— \gRBLVd 


\gRBLV u 


Mblr 
2 


2 v xr9bl 


~\ v xr9b 







M 2 


\gLv d 


-\gLv u 













— \9rblV(1 


\ghVd 







-\gRVd 










\9rbl v u 


-\9lV u 







\gRVu 










Mblr 
2 





-\gRVd 


\gRVu 


Mr 


-\ v xr9r 


\ v xr9R 




2 v xr9bl 











-\ v x R gR 





~^r 


V 


~\ v xr9bl 











\ v xr9R 


~Hr 






(42) 

The eigenvalues of this matrix are not completely arbitrary. Since U(l) b-l x U(1)r is broken 
in such a way as to produce correctly the SM group U(l)y in the limit of v << v R the matrix 
contains one state which corresponds to the MSSM bino, B, which is a superposition of Xbl 
and Xr. In addition the matrix contains an orthogonal state, which we will call B± in the 
following. 

For CMSSM like boundary conditions, M B l = M 2 = M R = Mx/ 2 , the bino is usually the 
lightest of the three gaugino like states, with the W being approximately twice as heavy. 
The B± is very often mixed with one of the right higgsinos, and, since Mbl at low energies 
is much smaller than vr this mixing is often important. In addition, there is the standard 
quasi-Dirac pair of "left" higgsinos, plus two more states which are mostly right higgsinos. 
Of the latter one is usually rather heavy, while the other can be light, if lir is small. 

In fig. the neutralino masses and TZ 2 ^ for CmBLR are plotted against vr for some 
arbitrary choice of other parameters. As discussed, there are in total seven eigenstates. Of 
special interest is the B±, so in the plot on the right we show the percentage of B± 
in the corresponding mass eigenstate. Here TZ^ = 1 means that the 2-th neutralino is a 
pure B±. As one can see in fig. |6]the masses and mixing of the three new states depend 
strongly on vr. For small vr all three states mix to each other. Increasing vr leads to 
a decoupling of the lighter higgsino-right from the B± which decreases in mass since hr 
becomes smaller for large vr while the masses of the two remaining states get large. Since 
the MSSM Neutralinos mix very little with the new states, there are four eigenvalues which 
show almost no dependence on the parameters vr and /ir. 

In fig. [6] the neutralino masses and 7Z 2 ±i are plotted against vr for the case of x^mBLR. 
In this calculation, fiR and vtia r can take fixed values while vr is varied freely. Two of 
the three new Neutralino states are a mixture of the higgsino-right and B± and therefore 
depend on vr. Since the lighter higgsino-right hardly mix to the Bj_ it has a constant mass 
at rnj lR ~ \fi R \ = 1700 GeV in this example. The lighter of the two new states that show 
dependence on vr is mostly a B±, whereas the one with larger mass is mostly a higgsino- 
right. The smaller vr the smaller the mixing between these two states and thus the larger 
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FIG. 5: Neutralino masses (left) and TZ 2 ±i (right) versus vr for otherwise fixed choice of parameters: 
m = 1000 GeV, M 1/2 = 1000 GeV, tan/3 = 10, A = -600, tan/3 R = 1.04. This plot uses the 
CmBLR version of the model. 
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FIG. 6: Neutralino masses (left) and K 2 ±i (right) versus Vr for otherwise fixed choice of parameters: 



m = 630 GeV, M 1/2 = 1000 GeV, tan/3 = 10, A = 0, tan/3 R 
itia r = 4800 GeV. This plot uses the XiJ m BLR version of the model. 



1.05, ii R 



■1700 GeV, 



the coupling of the mostly I?_|_-state to the MSSM particles. This will be important when 
we discuss LHC phenomenology in section IIVEI 

The dependence of the neutralino masses and TZ 2 ±i on p, R is shown in fig. [71 Since the 
higgsino-right and the B± mix, all three states show a dependence on /xr. The state which 
hardly mixes to the B± decreases in mass for small \hr\. So we can easily have a higgsino- 
right as LSP choosing /j,r close to zero. The state which is mostly the B± gets a smaller 
mass for large |/Ar|, while the one which is mostly a higgsino-right increases in mass. 



2. Sleptons and sneutrinos 



In models in which lepton number is broken, the scalar neutrinos split into a real and an 



imaginary part with slightly different masses 



76 



Since we assume that the smallness of 
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FIG. 7: Neutralino masses (left) and K 2 ±i (right) versus ^lr for otherwise fixed choice of parameters: 
m = 400 GeV, M 1/2 = 700 GeV, tan/3 = 10, A = 0, v R = 5000 GeV, tan/3 fi = 1.05, m AR = 
3000 GeV. 



neutrino masses is due to the smallness of the parameter fi s (and, therefore, _B Ms is supposed 
to be small too), this splitting between sneutrino mass eigenstates is too small to be of any 
relevance, except neutrino masses themselves. 

Neglecting fi S and B^ s the sneutrino mass matrix is given by 



M, 



2 



I ™Il,u ^Yjv u (A -cot(3fi) \v u v XR YlY s N 

^Y v v u {A - cot Pff) m\ R - ±Y s v XR (A - cot (3 Rf i R ) 

[ \v u v XR Y}Y v j-Y2v Xa (Ao- cot Prvlr) m| + '^Y}Y S ) 



(43) 

where 



v 2 



m LL,u = ™> 



I {(9bL + 9BLR ~ 9BLgRBL){v 2 XR ~ vl 2 ) + {g\ + g\ + QblQRBL) ~ V%)) 1 



v 2 + v 2 



RR,u — "V + 2 Y v Yv + 2 ^s^ s 



a Wbl + 9r + 9 2 blr + 9rbl ~ ^9bl9rbl - 2g R g B LR)(vl„ - v 2 



XR Xr> ' 

{9r + 9 2 RBL ~ 9BL9RBL ~ 9r9BLr) {v\ - 1%)) 1 ( 44 ) 

For charged sleptons one gets: 

M 2 = ( m l L J ^(A -tan^)\ 

' m 2 RR - ) 
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where 



m LL,f = m L + ^ Y l Y l ~ g (i&BL + 9 2 BLR ~ 9BLgRBL ~ 9r9BLR.) i. v \ R ~ ujj 

- (9l - 9bl9rbl - 9R9BLR){v 2 d - vl))l 

m RR,l = m E + Y YlY l + \ ~9 2 R + 9 BLR ~ 9rBl){^ r ~ «J H ) 

- (fl'l + fl'lsi + 9bl9rbl + 9r9blr) (v 2 d - vl)^j 1 (46) 

In fig. [Hlsneutrino and slepton masses are plotted against Vr, tan/3 R and li r . The figures 
on the left show a zoom into the region of the lightest states, whereas the figures on the right 
show a larger range of masses for a better understanding of the overall behavior. To see 
which particle is the LSP, while varying vr, tan /3r and lir, we included in all plots on the 
left the mass of the lightest neutralino state. This state is always a Bino, except for the plot 
against /Jr. Here the LSP becomes a higgsino-right for |/xr| < 250 GeV. The plots show that 
the masses depend strongly on the choice of vr and tan (3r. In the case of charged sleptons 
the dependence on vr and tan /3r comes only from additional D-terms at the tree level. This 
is different for sneutrinos. Here we can have an interplay between new D-terms and terms 
coming from the coupling Y s which both depend on vr and tan (3r. The additional D-terms 
force left sparticle to become light for tan /3r < 1 while for tan /3r > 1 right sparticle masses 
decrease. Up to vr — 6 TeV v\ is a right handed sneutrino and therefore the mass increase 
for increasing vr. For vr > 6 TeV the mass of v\ drops down again since here it is mainly 
a left handed sneutrino. Thus, increasing vr leads to a level-crossing in the mass spectrum 
of left and right handed sneutrinos. The same holds for the sleptons. In the plot against vr 
the mass of the right sneutrino decreases much faster for v r < 6 TeV than the mass of the 
right sleptons. This is due to the off-diagonal terms proportional to Y s , which contain also 
fiR in the sneutrino mass matrix. These terms mix the scalar component of S to vr. Thus 
for low values of vr in this example the LSP is neutral, which is allowed, whereas for larger 
values of vr (with left sleptons being light) there are parts of the parameter space, where 
the lightest slepton is charged, which is phenomenologically forbidden. Whether in the left 
sector charged or neutral states are lighter, depends heavily on the choice of parameters. 

Varying tan @r the right slepton masses decrease faster than the right sneutrino masses for 
tan/?R > 0.95 due the additional sneutrino mixing. Since the sneutrino and slepton masses 
depend strongly on the choice of vr, tan/3^ and lir one obtains limits on combination on 
these parameters. On the one hand, one has to avoid tachyonic states and on the other hand 
one has to take care not to get charged sleptons as LSP. The combination of both conditions 
forces us to choose tan (3r close to one and gives us an upper limit on vr and \[Ir\ as function 
of | tan /3r — 1 1 . 



19 




> 
o 



> 
o 





-1000 -750 -500 -250 
m [GeV] 



250 500 750 





1000 - 




800 - 


> 








O 


600 - 


5 


400 - 




200 - 



h 

h 



-1000 -750 -500 -250 
Hr [GeV] 



250 500 750 



FIG. 8: Lightest slepton (and neutralino) masses as function of vr, tan /3r and hr for a fixed 
but arbitrary choice of other parameters: tjiq = 220 GeV, M-y^ = 630 GeV, tan/3 = 10, Aq = 0, 
tan /3 R = 0.95, v R = 6000, fi R = -850 GeV, m AR = 2200 GeV, Y Sjii = 0.3. Plots on the left show 
a zoom into the light mass region, such that mass differences between the lightest sneutrino and 
the lightest charged slepton are resolved, figures to the right show the overall dependence, for a 
discussion see text. 



IV. CONSTRAINTS, SAMPLE SPECTRA AND DECAYS 

In this section we discuss several interesting phenomenological aspects which potentially 
allow the BLR model to be discriminated from the MSSM at the LHC and exemplify the most 
important features for a few study points. We include a discussion of the direct production 
of new states and characteristic changes in the cascade decays of supersymmetric particles. 

For brevity we will call these benchmark points BLRSP1- BLRSP5, the corresponding 



20 





BLRSP1 


BLRSP2 


BLRSP3 


BLRSP4 


BLRSP5 


cMSSM 


m [GeV] 


470 


1000 


120 


165 


500 


M 1/2 [GeV] 


700 


1000 


780 


700 


850 


tan /3 


20 


10 


10 


10 


10 


A 





-3000 


-300 





-600 


Extended gauge sector 


vr [GeV] 


4700 


6000 


6000 


5400 


5000 


tan/3/j 


1.05 


1.025 


0.85 


1.06 


1.023 


/x B [GeV] 


-1650 


-780 


-1270 


260 


(-905) 


m AR [GeV] 


4800 


7600 


800 


2350 


(1482) 


Yukawas 




0.04 


0.1 


0.1 


0.1 


0.1 




0.04 


0.1 


0.1 


0.1 


0.1 


^,33 


0.04 


0.1 


0.1 


0.1 


0.1 


n,H 


0.04 


0.042 


0.3 


0.3 


0.3 


V s , 22 


0.05 


0.042 


0.3 


0.3 


0.3 


^s,33 


0.05 


0.042 


0.3 


0.3 


0.3 



TABLE II: Parameters of the various study points. In BLRSP1-BLRSP4 /xr and rriA R are input 
whereas in BLRSP5 the constrained version of the model has been used and, thus, these two 
parameters are output. For a discussion of these points see text. 



input parameters are listed in table [UJ All of these points have been calculated with the 
XijHiBLR version of the model. However, note that for BLRSP5 the input is chosen to be 
consistent with the CmBLR variant. 

A few comments on the input parameters and the resulting mass spectra are in order, 
before we discuss the phenomenology in detail. As shown below, the bounds on rare lepton 
flavour violating decays require Y u and Y s to be essentially flavour- diagonal, unless these 
couplings are very small. Therefore we have chosen Y u and Y s diagonal as starting point 
implying that all points satisfy trivially the LFV constraints. A correct explanation for the 
neutrino angles then requires flavour violating entries in the parameter /is, which we do not 
give in table [Til since they are irrelevant for collider phenomenology. 

The input values of table [TTJ lead to the mass spectrum shown in table IIHI We give 
the masses and in brackets the particle character. In case of mixed states the two largest 
components, for example (W, Hl), are given where the first entry accounts for the larger 
contribution. If the ordering in the composition changes like in the case of m^ 5 6 we use 
squared brackets. Therefore we have (cl,ul) for m u 5 and (ul,cl) for m Sg . In all cases 
input parameters have been chosen such, that the squark and gluino masses are outside the 



region currently excluded by pure CMSSM searches at ATLAS [44j and CMS [45|. Since 
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(a) we expect the missing momentum signal to be smaller in these points than in a true 
CMSSM spectrum and (b) our squark spectra are less degenerate than the CMSSM case, 
we believe this is a conservative choice. Two of the points have a sneutrino LSP (BLRSP1 
and BLRSP3), while three points have a neutralino LSP (for BLRSP2 and BLRSP5 mostly 
a bino, for BLRSP4 a state which is mostly a Hr). 

Note that the ordering of sfermion mass eigenstates does in many cases not follow the 
standard CMSSM patterns: < mp, R ~ < mp, L ~ mg L < mf 2 and m~ tl < rric R — 
rriQ R < m>c L — < mr (similar for sdowns). These patterns are distorted in the study 
points due to the unconventional D-terms of the model and this feature gets enhanced for 
larger | tan/3_R — 1| and/or larger values of vr. We note also that for sneutrinos and charged 
sleptons many states are quite degenerate. For example fiR and cr have practically the same 
mass in all points. While these degeneracies are always true in CMSSM spectra, in our case 
this is not necessarily so, but simply reflects the fact that both Y v and Y s have been chosen 
generation independent in all points, except BLRSP1. As this point shows, even a rather 
moderate generation dependent value of Y s can lead to large mass splittings in the sneutrino 
sector. A generation dependent value of Y v would not only split sneutrino masses but also 
charged slepton masses. 



A. Lepton decays and LFV 



To explain the measured neutrino mixing angles by the mass matrix given in eq. ( 138]) 
the Yukawa couplings Y v , Y s and/or the the bilinear term fi$ have to contain off-diagonal 
elements. In case that the neutrino mixing is explained by the form of Y v or Y s also lepton 
flavor violation in the charged lepton sector will be induced. On the one hand the contri- 
butions of Y v to the RGE evaluation branching ratios of Y e and to the soft-breaking terms 
in the lepton sector open decay channels like li — > Ijj and /j — > 3lj similarly to high-scale 
seesaw type I — III 77H79I] . On the other hand, in inverse seesaw the entries of Y v can be po- 
tential large and the dominant contributions to LFV decays can come from diagrams which 
are proportional to (Y u Yj)ij. For a long time it has been assumed that the most stringent 
bounds on (Y u Yj)ij come from the radiative decay fi — > e'-f while the photonic contributions 
to n — y 3e are always smaller and therefore Br(/i — > e'-f) > Br(/i — > 3e) must hold. How- 
ever, recently it has been pointed out that in presence of new Yukawa couplings like Y v the 



Z-penguin contributions to Br(yU — > 3e) can dominate 80] . These are less suppressed then 
the photonic contributions by a factor 



MSUSY 

M z 



As consequence, the experimental limits 
on the decay into two charged electrons can be much more constraining than the radiative 
decay in case of a heavy SUSY spectrum. We can see this by parametrizing the neutrino 
Yukawa coupling as 



(47) 
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1247.0 (tR,t L ) 


1111.3 (tfl.tx) 


1316.0 (tR,f L ) 




[GeV] 


1392.1 (i L) t fl ) 


1851.9 (f L ,tR) 


1526.9 (t L ,t R ) 


1361.4 (f £j f fl ) 


1643.2 (f L ,tH) 


m «3,4 


[GeV] 


1456.0 {cr)/(ur) 


2154.7 (cr)/(ur) 


1565.9 (c R )/(u R ) 


1392.4 (cr)/(ur) 


1728.0 (cr)/(ub) 


m U5,6 


[GeV] 


1509.0 [c L ,fi L ] 


2227.3 [c L ,u L ] 


1634.0 [c £ ,u £ ] 


1448.8 [c Lj u L ] 


1795.8 [c L ,u L ] 


m d, 


[GeV] 


1359.2 (b L ~b R ) 


1819.2 (6 L ) 


1409.8 (b R ~b L ) 


1326.3 (b L ) 


1611.8 (6 X ) 


m d 2 


[GeV] 


1464.0 (6r,6 l ) 


2148.1 (6b) 


1462.3 (sr) 


1420.1 (6r) 


1724.5 (6«) 


m d 3 


[GeV] 


1489.8 (S fl ) 


2175.9 (Sr) 


1462.3 (da) 


1426.2 (S R ) 


1734.8 (S R ) 




[GeV] 


1489.8 (d R ) 


2175.9 (4) 


1496.2 (b L ~b R ) 


1426.2 (d fl ) 


1734.8 (Jr) 


mi 

0.5,6 


[GeV] 


1509.0 [S L ,d L ] 


2228.9 [s L ,d £ ] 


1635.9 [s £ ,d L ] 


1450.9 [a L ,rJi] 


1795.8 [s L ,d L ] 


Higgs (l-loop/2-loop) 




[GeV] 


59.1/59.6 


119.2/125.4 


92.7/93.1 


100.8/102.6 


18.8/18.8 


m h2 [GeV] 


119.0/124.1 


139.7/140.4 


114.5/120.1 


121.0/124.8 


115.7/121.8 




0.05/0.04 


0.90/0.83 


0.07/0.04 


0.33/0.22 


0.001/0.001 




0.95/0.96 


0.10/0.17 


0.93/0.96 


0.67/0.78 


0.999/0.999 


Neutralinos 




[GeV] 


282.2 (B) 


416.7 (B) 


312.9 (B) 


258.5 (h R ) 


346.6 (B) 




[GeV] 


552.3 (W,h L ) 


780.0 (h R ) 


615.3 (W", fo) 


279.7 (B) 


679.5 (W,h L ) 




[GeV] 


828.9 (fo) 


817.5 (W") 


1086.6 (ft £ ) 


549.0 (W'.fcr) 


902.7 (h R ) 




[GeV] 


838.9 (h L ,W) 


1865.5 (hi) 


1092.8 (h L ,W) 


844.9 (h L ) 


1133.1 (h L ) 


m x o 


[GeV] 


1230.4 (B^r) 


1865.7 (h L ) 


1232.2 (h R ,B ± ) 


856.8 (h £ , W") 


1139.4 (h i; W) 


m x£ 


[GeV] 


1650.9 (/ir) 


2017.6 (B ± ,4) 


1811.3 (B ± ,4) 


1639.0 (Bj_, 6r) 


1489.8 (B ± ,h R ) 


m x a 7 


[GeV] 


2608.3 (h R ,B±) 


2392.3 (h R , B±) 


2741.4 (4,B ± ) 


2174.6 (4,B ± ) 


2056.5 (h R ,B ± ) 



TABLE III: Spectra of our study points, for discussion see text. (v R ) is a nearly maximal mixture 
of the right sneutrinos and the 5-fields. 
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FIG. 9: Branching ratios of lepton flavour violating processes as a function of / for tan j3 = 10, 
A = 0, v R = 5000 GeV, tan fa = 1.05, fi R = -500 GeV, m AR = 1000 GeV and three (m ,M 1/2 ) 
combinations. 



with 

a = (Am!/ Ami)* ~ 0.4 and 6 = 0.23. (48) 

Here, Am^ and Am\ are the mass differences measured in solar and athmospheric neutrino 
oscillations. The value of b accommodates for sin 2 # 13 = 0.026. 

Br(/i — > ej) and Br(/i — > 3e) as function of / are depicted in fig. El and we show also the 



most recent, experimental limits of 51 



Br(/x -> erj) < 2.4 • 10" 12 , Br(/i -> 3e) < 1.0 • 10" 12 (49) 

For the plot we chose three different points for (mo, Mi/ 2 ). Since Br(/i — >■ 3e) hardly 
depends on SUSY masses for msuSY 3> all three lines lie very close together in contrast 
to Bt(h — > ej). For light SUSY spectra Br(/i — > e^y) is dominant whereas the heavier the 
SUSY particles the more important gets Br(yU — > 3e). As shown in BLRSP1 and BLRSP2 
we can have light right-handed neutrinos such that contributions from the H^-graph to LFV 
can not be neglected anymore. In fig. [10] Br(yU — > c-f) and Br(/i — > 3e) are plotted as a 
function of m VR . For masses below 300 GeV contributions from right-handed neutrinos start 
to dominate. The minimum in Br(/i — > 3e) comes from a cancelation between the right- 
handed neutrino and the corresponding SUSY graph. In the limit of large m VR Br(/i — > cy) 
and Br(/i — > 3e) converge to the value coming from SUSY contributions. 

However, it is possible to circumvent these bounds by assuming that Y v and Y s are 
diagonal and the entire neutrino mixing is explained by fi$- Of course, this will not only 
reduce Br(/x — > 3e) but also Br(yU — > ej). 



B. Higgs physics, direct production 

In all study points of table [TT] there is one Higgs boson with mass between 120 and 
125 GeV. In addition, there is a second state with masses varying between 19 and 140 GeV. 
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10" 9 

10-10 
T lo- 11 

PQ io- 12 
lO" 13 

500 1000 1500 2000 2500 3000 3500 4000 

m VR 

FIG. 10: Branching ratios of lepton flavour violating processes as a function of m VR for niQ = 
800 GeV, M 1/2 = 1200 GeV, tan/3 = 10, A = 0,v R = 10 TeV, tan(3 R = 1.05, (jl r = -500 GeV, 
rriA R = 1000 GeV. The dashed red line is the upper limit for /i — > e'y and the dashed blue line for 
r — > hj. 




In BLRSP1, BLRSP3 and BLRSP5 the mass eigenstate h 2 is SM-like, with R\ 2 > 0.9. In 
BLRSP2 it is hi, which has a large content of Hd and H u and BLRSP4 is a case where 
hi and hi have large mixing. Since we have often a mass eigenstate below the LEP limit 
of 115 GeV for a standard model Higgs boson, we have checked the consistency of these 



eigenstates with data using HiggsBounds 3.4.0beta |82l . [83j. All points are allowed by 
accelerator constraints, but sometimes very close to existing bounds, especially BLRSP4 
and also BLRSP2. As an indication for the theoretical uncertainties in the mass calculation 
we give the masses using the complete 1-loop formulas and the ones adding the dominant 



2-loop corrections to the MSSM sector 8J-|89 



In BLRSP1 and BLRSP5 hi is so light that the decay h 2 — > hihi is kinematically allowed. 
However, the mixing between both sectors is so small that for BLRSP1 the corresponding 
branching ratio is about 1 per-cent whereas for BLRSP5 it is a few per-mile. The smallness 
of this decay is a direct consequence of the bounds imposed by LEP and the decay h 2 — > hihi 
can never be dominant in the BLR model. The h 2 can also decay into a combination of 
heavy and light neutrinos with a branching ratio of a few per-cent, as for example in case 
of BLRSP1 leading to the final states 

h 2 -> UiU k -> v^W* (50) 
h 2 -> ViV k ->■ ViVjZ (51) 

with i,j = 1, 2, 3 and k — 4, ... , 9. These final states can also be obtained via intermediate 
states containing an off-shell vector boson, e.g. WW* and ZZ*. However, their existence 
implies that ratio of quark versus lepton final states will not correspond to the branching 
ratios of the vector bosons. Note, however, that for hadronic W-boson decays the invariant 
mass of jj+lepton system would show a peak at the heavy neutrino mass, which allows to 
identify this signature, in principle. Apart from these decays, the h 2 can also decay to two 
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scalar neutrinos and, if kinematically allowed, this decay can become dominant, leading to 
a (nearly) invisible Higgs boson. 



The hint for a 125 GeV Higgs boson [3, [8|, see also |44| and |45j, indicates a slightly larger 
than expected branching ratio into the two-photon final state. In 18j it was shown that 
the NMSSM can, in principle, explain such an enhanced di-photon rate, due to a possible 
mixing of the singlet and the Higgs, which reduces the coupling of the Higgs to bottom 
quarks, thus reducing the total width, without affecting the production cross section. In the 
case of the BLR model, such a construction is not possible, since our singlets are charged 
under U(1)r and the mixing between SM and BLR sectors is controlled by tan/3^. Since 
we have to choose tan/3j{ close to one, the singlets mix to the up and down components of 
the Higgs equally. Therefore a reduction of h — > bb causes a reduction of the coupling for 
gluon fusion as well. Thus, a sizeable enhancement of Br(/i — > 77) by reducing simply the 
total width is not possible in the BLR model. Currently the discrepancy of the data with 
expectations is only at the level of about 1 a c.l. However, should future data show indeed 
an enhanced rate for the 77 final state, this would be hard to explain in the BLR model. 

In the four points (BLRSP1, BLRSP3-BLRSP5) hi has approximately the same branching 
ratios for the decays into SM-fermions as a SM Higgs boson of the same mass. However, 
the corresponding widths are suppressed by the mixing with the usual MSSM sector which 
reduces the width by a factor between 10 2 and 10 4 . At the LHC the main production of 
this particle is via SUSY cascade decays, e.g. it appears in the decays of £4,5,6 (BLRSP1, 
BLRSP3), X3 (BLRSP4) or in the decays of the heavy neutrinos which are produced via the 
Z' (BLRSP1, BLRSP4) as discussed in section [TV] However, in case of BLRSP5 LHC will 
miss h\ as it only appears in the decays of the heavy Higgs bosons which have masses in the 
TeV range. 

Study point BLRSP2 differs from the others as here hi is the MSSM-like Higgs boson 
and h 2 has a mass of 140 GeV which could explain the slight excess in this region observed 



by ATLAS and CMS in the early data 9C| |9l| . In this region of the parameter space the 



Higgs at 125 GeV is made as in the MSSM, implying a rather heavy SUSY spectrum. This 
is due to the fact that a 140 GeV Higgs with reduced couplings can only be the hsLR, i-e. 
this points exist to the right of the level-crossing region shown in fig. [2j Due to the choice 
of a rather small Y s in this point the heavy neutrinos masses are below the mass of h 2 . This 
leads to non-standard decays into the heavy neutrinos which dominantly decay to a lepton 
and a W-boson. 



C. Z' physics 



As already mentioned in sect. IIII A\ our Z' corresponds essentially to the Z x in the nota- 



tion of 



75 



In previous studies usually two assumptions have been made in the construction 
of mass bounds: (i) the Z' decays only into the known SM particles 92[ and (ii) the effects of 
gauge kinetic mixing are neglected. Both assumptions are not truly valid in the BLR model. 
As shown in table |IV] we find in all our points that the heavy neutrinos appear as final states 



26 



TABLE IV: Branching ratios of the dominant Z' decay modes. Here we have summed over the 
generations in case of the charged fermions and sfermions. For the neutrinos we have splitted this 
sum into a sum over the light (heavy) states denoted by i>\ iy^)- 

final state BLRSP1 BLRSP2 BLRSP3 BLRSP4 BLRSP5 



BR(dd) 


0.31 


0.35 


0.35 


0.37 


0.43 


BR(uu) 


0.06 


0.07 


0.07 


0.07 


0.08 


BR{ll) 


0.12 


0.14 


0.14 


0.14 


0.16 


BR{u m ) 


0.10 


0.11 


0.12 


0.12 


0.12 


BR(v h v h ) 


0.27 


0.30 


0.13 


0.11 


0.13 


BRivv) 


0.05 




0.05 


0.03 




BR(!l) 






0.05 


0.03 




BR(xtx2) 








0.02 




BR(xix° 5 ) 








0.02 







CL 


CR 


d 


-i (-3g L Z 13 + g BL Z 23 + g BLR Z 33 ) 


-£ {(sbl - 3g RBL )Z 23 + (g B LR - 3g R )Z 33 ) 


u 


-i {3g L Z 13 + g BL Z 23 + g B L R Z 33 ) 


-| {(9BL + 3g RBL )Z 23 + (g BLR + 3g R )Z 33 ) 


1 


| (g L Z 13 + g BL Z 23 + 9B lrZ 33 ) 


{ {(9bl + 9rbl)Z 23 + (g BLR + g R )Z 33 ) 


V 


J [ELi Zi 3+x '*Z*+*{(-g BL + g RBL )Z 23 
+{9R- 9blr)Z 33 ) + 
ELi Zt' ,m Z? (9blZ 23 - g L Z^ 3 + g BLR Z 33 ) ' 


-|[Lx=l Z ^ z » {{-9BL+ 9rbl)Z za 
+{9R- 9blr)Z 33 ) + 
Ex=l Zi x '*Zt (g BL Z 23 - g L Z^ 3 + g BLR Z 33 ) ' 



TABLE V: Coefficients c L and c R for the coupling between Z R and two leptons or quarks. Here, 
Z is the rotation matrix diagonalizing the neutral gauge boson mass matrix and Z v is the neutrino 
mixing matrix. 



beside the SM-fermions. Moreover, in all but BLRSP5 also supersymmetric particles appear 
as decay products, in particular sneutrinos and sleptons. On the other hand, gauge kinetic 
effects are in this model less important and were only important if one could measure the 
branching with a precision of 1 per-cent or better. 
The Z' couples to leptons and quarks as follows 

Z'Jr{c{P L + c R P R )f (52) 

The different coefficients are given in table |V] 

Note, that in the couplings to the u-quarks a partial cancellation occurs in contrast to 
the ones to <i-quarks, which get enhanced. Moreover, the same feature appears in the vertex 
q-q-B± which leads to some interesting consequences discussed in section IIVEI 

We find that the decays into the heavy neutrino states are always possible and have a 
sizable branching ratio provided Tr(|y s |) < 1. In table [IV] we summarize the most important 
final states of the Z 1 for the different scenarios. As can be seen the heavy neutrino final 
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b L , 1 , , , , 1 , , , , 1 , , , , L_ 

1500 1600 1700 1800 

m z , [GcV] 



FIG. 11: Cross section of pp — > Z' — > p. + p>~ near the Z' peak as function of mz> taking into 
account a K-factor of 1.3 [98(. For the black, dotted line the Z 1 width has been calculated allowing 
only SM final states, while the blue solid includes also right-handed neutrinos and SUSY states. 
The red line shows the ATLAS exclusion limit 



74i | . We have used as input BLRSP1 and varied 



v R = [4.1,5.1] TeV. 



states have always a sizable branching ratios with up to about 30 per-cent when summing 
over the generations. But even for rather heavy neutrinos as in BLRSP5 own finds for this 
channel a 15 per-cent branching ratio. In several cases also channels into SUSY particles are 
open, in particular in scenarios with sneutrino LSPs. In case of supersymmetric particles 
the final states containing sleptons or sneutrinos have the largest branching ratios. Channels 
into neutralinos or charginos are suppressed. They proceed either via the mixing with the 
Z which is rather small or via the projection of the higgsino-right onto the corresponding 
neutralino state. 

The appearance of additional final states leads to a reduction of the event numbers in the 
most sensitive search channels, i.e. reducing cross section times branching ratio, and, thus, 



the bounds obtained by the LHC collaborations [74|, |93|, |94j are less constraining in the BLR 
model. This is depicted in fig. [TT1 where we show the production cross section a(pp — > 
arround the Z' resonance. 7 In case that the width of the Z' is calculated using only SM 
final states the cross section is increased roughly by a factor 1.6 in comparison to the case 
where also right handed neutrinos and SUSY particle contribute to the width of Z'. With 
this choice of parameters, the main effect is due to R-neutrinos. We attribute the remaining 
difference to the official ATLAS result to slightly different values in the couplings and slig htly 



different branching ratios of the final states. Our results agree also with the ones of ref. [75 
We conclude that, although in our benchmark points we take always vnz> > 1.8 TeV, a 
significantly lower mass is possible consistent with data. 



7 For the calculation of the cross section we used WHIZARD 
SUSY Toolbox b" 
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and implemented the model using the 
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D. Heavy neutrinos 



As discussed above, the heavy neutrino states can be produced via the Z' with a con- 
siderable branching ratio of about 30 per-cent when summing over all heavy neutrinos. 
Moreover, see below, they can also be produced in the cascade decays of supersymmetric 
particles. These heavy neutrinos mix with the light neutrino states implying a reduction 
of the couplings of the light neutrinos to the Z-boson and, thus, also a reduction of the 
invisible width of the Z-boson. Taking the data from ref. 5l| this can be translated into 
the following condition on the 3x3 sub-block U^, i,j < 3, of the neutrino mixing matrix: 



(53) 



3 




2 










< 0.009 


y=ii*<j 


k=l 







at the 3-cr level. We have checked that all our benchmark points fullfill this condition. 
The main decay modes of the heavy neutrinos are 8 



Zvi 



(54) 
(55) 
(56) 



where j > 4, % < 3, k — 1, 2 and I = e,fi,r, provided they are kinemtically allowed. 
If there is no kinematical suppression we find in general the branching ratios scale like 
BR(uj ->■ W ± l T ) : BR(vj ->■ Zv^) : BR(uj -> h k Vi) ~ 0.5 : 0.25 : 0.25 where we have 
summed over the light Higgs bosons, the light neutrinos and leptons, respectively. We stress 
that these states are quasi-Dirac neutrinos and, thus, for six heavy neutrinos at LHC the 
existence of up to three new particles could be established. Note, that the final states 
containing a PU-boson allow for a direct mass measurement. 

Beside the above decay modes, also decays into SUSY particle are possible if kinematics 
allow for it. For example we find that for BLRSP4 the decay into ^1,2,3X1 are possible and 
have branching ratios of about 3 per-cent. In scenarios like BLRSP3, BLRSP4 and BLRSP5 
the main production of the heavy neutrinos is via the Z' and, thus, a high luminosity will 
be required to observe such final states. 



E. SUSY cascade decays 

In this section we point out several features which distinguish the BLR model from the 
usual MSSM. For the sake of preparing the ground, let us first summarize the main features 
of the MSSM relevant for the LHC, focusing for the time being on scenarios where the gluino 
is heavier than the squarks: (i) The gluino decays dominantly into squarks and quarks, (ii) 



For related discussions see e.g. [99l4l0l| and references therein 
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L-squarks and L-sleptons decay dominantly into the chargino and the neutralino which are 
mainly S , f/(2)i-gauginos. Apart from kinematical effects the branching ratio for decays into 
the charged wino divided by the branching ratio into the neutral wino is about 2:1. (iii) 
R-squarks and R-sleptons decay dominantly into the bino-like neutralino with a branching 
ratio often quite close to 100 per-cent. (iv) In case of third generation sfermions also decays 
into higgsinos are important. 

In the BLR model one has two main new features: (i) there are additional neutralinos 
and (ii) the sneutrino sector is enlarged as well. The latter implies that sneutrino LSPs 
are possi ble consistent with all astrophysical constraints and direct dark matter searches 
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1081 ] . This feature is for example realised in study points BLRSP1 and BLRSP3 . 
Let us start the discussion with BLRSP1. In this point the four lighter neutralinos are 
the usual MSSM neutralinos with the standard hierarchy. The fifth state corresponds to 
the additional £/(l)-gaugino, which we call B±, whereas the two additional states are the 
additional higgsinos. Note that the lighest neutralino is not stable anymore but decays 
into final states containing all nine neutrinos as well as the three lightest sneutrinos. Of 
the latter ones the second lightest is so long lived that it will lead to a displaced vertex 
in a typical collider detector. The third sneutrino decays dominatly via three-body decays 
into l + l~b>i and v^i^i with % = 1,2 and k, I = 1,2,3. As discussed in section TlV Dl the 
heavy neutrinos decay dominantly into W^-bosons and charged leptons, thus the decays of 
the lightest neutralino are not invisible. 

Bj_ appears for example in the decays of and sr with branching ratios BR(gx°) ~ 0.8 
and BR(gXg) ~ 0.2. For completeness we remark that the decays of ur and cr into x® is 
suppressed as the corresponding coupling is supressed as are the couplings of Z' to u-type 
quarks in this model, xt decays dominantly into sleptons and sneutrinos. Combining all 
the above together one gets a much richer structure for the decays of the R-squarks, e.g. the 
following decay chains: 

Qr -> QXi ~> QVkVi ->■ qVjZvx (57) 
Qr qXi -> qv k vi ->■ ql ± W T u 1 (58) 

q R qx°i ->• q»kh ->• q^wn'+i'-u, (59) 
qR -> qxl -»• q^if -> qi^xl ->• qPFwi qi ± i T i' ± w T D 1 (60) 

with k G {4, 5, 6, 7, 8, 9} and j G {1, 2, 3}. Of course, several other combinations are possible 
as well. 

From equations ( 1571) to ( 160]) one sees immediatly that the standard signature of R-squarks, 
namely jet and missing energy, is only realized in a few cases in this study point, e.g. if in 
eq. ( 15"7|) the Z decays into neutrinos. Interestingly, the chain via x\ m to sleptons leads to 
a characteristic edge in the invariant mass of the lepton which can be used to determine 



9 We note for completeness, that the relic abundance is actually somewhat too large in this point but can 
easily be adjusted by changing for example in BLRSP1 tan/3^ from 1.05 to 1.0475 without changing the 
collider features of BLRSP1. 
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the corresponding masses once combined with information from other decay chains. Also 
in the study points BLRSP2 and BLRSP5 cIr and §r decay into heavy neutralinos, which 
contain sizable content of the extra U(l) gaugino, with a sizable branching ratio. However, 
there the situation is somewhat less involved as in these study points the lightest bino-like 
neutralino is the LSP. 

Another interesting feature is, that xt decays also into the heavier sneutrinos which 
themselves decay into the LSP plus h\. Similarely hi can be produced in the decays of the 
heavy neutrinos implying that this state can be produced with sizable rate in SUSY cascade 
decays. However, as the corresponding final states are quite complicated a dedicated Monte 
Carlo study will be necessary to decide if this is indeed a discovery channel for hi. 

From the point of view of SUSY cascade decays BLRSP2 looks essentially like a standard 
MSSM point. Inspection of the spectrum shows that x 2 i s essentially a higgsino correspond- 
ing to the extended U(l) sector but it shows hardly up in the cascade decays. Its main 
production channel is via an s-channel Z' but even in this case the corresponding cross 
section is so low that it will not be dedected at the LHC even with an integrated luminosity 
of 300 fb _1 . Another interesting feature shows up in the decays of X3 which is mainly the 
neutral wino and gets copiously produced in the decays of the L-squarks: it decays with 
about 77 (15) per-cent into h\ (h 2 ), implying that the cascade decays are an important 
source of Higgs bosons. 

In case of BLRSP3 one has sneutrino LSPs like in BLRSP1 but with a different hierar- 
chy in the spectrum, as the three lightest sleptons are lighter then the lighest neutralino. 
Therefore the Xi has also sizable decay rates into charged sleptons which sum to about 30 
per-cent. The sleptons decay then further into W~vi t2 ,3 and v 2 $ via 3-body decays into 
//-pairs. The latter, however, are rather soft due to the small mass difference. In addition 
we have the decay channel into a light neutrino and one of the heavier sneutrinos which 
themselves decay into a lighter sneutrino and either one of the Higgs boson or the Z-boson. 
Putting again all these decays together one obtains for the Xi decays 



X? 


-> 


i ± w T i>i 


(61) 


x? 


-> 


i±F -»• i ± w T i> 2:6 -»• i ± w T ffh 


(62) 


X? 


-> 


Vjf>2,3 -> ^1,2,3/M 


(63) 


X? 


-> 


Vfii 


(64) 


X? 


-> 


VjV k -> h>jhi )2 vi 


(65) 


*S 


-> 


VjV k ->> Vjhi^ffvi 


(66) 



with j = 1,2,3 and k = 4,5,6. This implies that the decays of the R-squarks show again 
a more complicated structure compared to the usual CMSSM expectations. Channels ( 165]) 
and (166]) give hi in about 15 per-cent of the final states of Xv Moreover, x 2 an d Xi decay 
dominantly into sleptons and sneutrinos. Here a new feature is found for xf, as also the 
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following chains 



(67) 
(68) 



gives rise to sharp edge structures. However, as the main final states of Z and h\$ are two 
jets, the feasability still needs to be investigated. 

In BLRSP4 we have chosen fiR = 260 GeV in order to construct an LSP which is essen- 
tially a Hr. Here, the -R-sleptons are lighter than Xz-* which is essentially bino-like in this 
point, giving rise to the following decay chain of the down-type i?-squarks 



Nearly all cascade decays end in a xl or one °f the lighter sleptons. Due to the fact, that in 
this particular case the additional sneutrino states are hardly produced, it might be difficult 
to disthinguish it from the NMSSM, at least as long as the Z' is not discovered. The heavier 
L-sleptons do not show up in the cascade decays of squarks and gluinos but can be produced 
via the Z' as discussed in section IIV1 

BLRSP5 is similar to BLRSP1 but compatible with pure GUT conditions, e.g. \Ir and 
rriA R are not input in this case put derived quantities. To fullfill the tadpole equations we 
have to choose Y s = 0.3 and tan /3r = 1.03 if we want a relatively low mo = 500 GeV 
while Mi/2 — 850 GeV. The choice of Y s leads automatically to large masses for the heavy 
neutrinos such that the lightest Higgs can not decay into those states. As in BLRSP1 the 
down-type i?-squarks decay not only into Xi but also into xt with a branching ratio of about 
13 per-cent. For completeness, we note that here Xz — h>R- However this state gets hardly 
produced in any of the SUSY decays or via the Z'. Therefore, it is likely that LHC will miss 
it and also at a linear collider such as ILC or CLIC it will be difficult to study, due to the 
small production cross section. 

V. CONCLUSIONS 

We have studied the minimal supersymmetric U(l) b-l x U (1)r extension of the standard 
model. The model is minimal in the sense that the extended gauge symmetry is broken with 
the minimal number of Higgs fields. In the matter sector the model contains (three copies 
of) a superfield z> c , to cancel anomalies. Adding three singlet superfields S allows to generate 
small neutrino masses with an inverse seesaw mechanism. 

The phenomenology of the model differs from the MSSM in a number of interesting 
aspects. We have foccused on the Higgs phenomenology and discussed changes in SUSY 
spectra and decays with respect to the MSSM. The model is less constrained then the 
CMSSM from the possible measurement of a Higgs with a mass of the order of 125 GeV. If 
the hints found in LHC data 0, 0] is indeed correct our model predicts two relatively light 
states should exist, with the second h° corresponding (mostly) to the lightest of the "right" 
Higgses, added to break the extended gauge group. 



d R -» dx° 2 -> dlH* -> d^l^Xi 



(69) 
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It is interesting, as we have discussed, that very often a right sneutrino is found to be the 
LSP. This will affect all constraints on CMSSM parameter space derived from constraints on 
the dark matter abundance. In fact, if the right sneutrino is indeed the LSP in our model, 
no constraint on any CMSSM parameters can be derived from DM constraints. 

The model has new D-terms in all scalar mass matrices, which can lead to sizeable 
changes in the SUSY spectra, of potential phenomeno logical interest. We have discussed a 
few benchmark points, covering a number of features which could allow to distinguish the 
model from the CMSSM. Obviously this includes the discovery of a Z' at the LHC where we 
have shown that the current bounds from LHC data depend on the details of the particle 
spectrum. Also the cascade decays of supersymmetric particles can be significantly more 
involved than in the usual CMSSM as the additional neutralinos, neutrinos and sneutrinos 
lead to enhancement of the multiplicities in the final states. This implies that the existing 
limits on the CMSSM parameter space get modified as standard final states have reduced 
branching ratios and at the same time additional final states are present. In case that the 
mBLR model is indeed realized these new cascade decays will offer additional kinematical 
information on the particle spectrum. 
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Appendix A: Appendix 
1. Mass matrices 

Here we list the tree-level mass matrices of the model not given in the main text. 



• Mass matrix for Down-Squarks, Basis: ( di, df> 





(Al) 
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m LL = m l + Y Y d Yd ~ ^ {{dBL + 9 2 BLR ~ 9BLR9R ~ 9bl9RBl) [f XR ~ u| H ) 

+ (3^1 + 9bl9rbl + 9blr9r) (yj - v 2 ) ) 1 ( A2) 



2 2 



24 ((^II + 9 blr ~ ^(9blr9r + 9bl9rbl) + 3(^ + 

+ [9BL9RBL + S-BLi^il ~ 3(^| + g-| si ) J - ^) J 1 (A3) 



Mass matrix for Up-Squarks, Basis: (ul,ur) 

m « = 1 / V * _j_ T ^ 2 I 

\72{~ VdYu ^ + v u T u) m RR 

m LL = m 2 q + Y Y ^ Yu ~ 2I ((^ L + ^« ~ ~ 9bl9rbl) (v 2 Xr - v\ R 

+ ( 3 #! - 9bl9rbl ~ 9blr9r) (y 2 u ~ v d)) 1 ( A5 ) 



2 

24 [\9bl + 9 blr + ^{9blr9r + 9bl9rbl) ~ 3(#| + 9rbl)) 

+ \9bl9rbl + 9blr9r + 3 (^ + [ v d - V u) ) 1 ( A6 ) 



■v 2 - 

XR Xr 



Mass of the Charged Higgs boson: One obtains the same expression as in the 
MSSM: 

m 2 H+ = (tan (3 + cot (3) + m\y (A7) 
Mass matrix for Charginos, Basis: {w~ ,H~^ , (w + ,H^j 



^9LV d fi 



2. Calculation of the mass spectrum 



We are going to present now the basic steps to calculate the mass spectrum. As starting 
point we use electroweak precision data to get the gauge and Yukawa couplings: the SM- 
like Yukawa couplings are calculated from the fermion masses and the one-loop relations 



of ref. 68| which have been adjusted to our model. Similarly, also the standard model 



gauge couplings are calculated by the same procedure presented in ref. 68], but again 



including all new contributions of the mode under consideration. Since the entire RGE 
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running is performed in the basis SU(3)c X SU(2)l x U(1)r x the value of the 

GUT normalized g^i and gu are matched to the GUT normalized hypercharge coupling gy 
by 

9r =cry gY , (A9) 



9bl 



^9blr9rbl9r ~ VQgRBLgy + V (3^blh _ <1 ^9blr9r + 2g R )(h(g R + g-| BX - 3^)^ 



%1 - 3 9y 

(A10) 



This is nothing else then an inversion of the well known relation between the gauge couplings 
for U(1)r x U(1)b-l U{1)y including the off-diagonal gauge couplings given in eq. flAHI) . 
We are using the 50(10) GUT normalization of for t/(l)y and for U(1)b~l- To 
get the correct values of cry as well as gRBL and gsLR an iterative procedure is used: cry is 
calculated as ratio of the gy and qbl when running down from the GUT scale and applying 



5(9bl9r — 9blr9rbl) 2 / A11 x 

9y = \l 7, t, 7, t= (All) 

%Il + 9blr) + 2{g% + g 2 RBL ) - 2^Q{g R g BLR + 9bl9rbl) 

When the gauge and Yukawa couplings are derived, the RGEs are then evaluated up 
to the GUT scale where the corresponding boundary conditions of eqs. (jSJ), (EJ) and ( I12j) 
are applied. Afterwards a RGE running of the full set of parameters to the SUSY scale 
is perfor med. We use always 2-loop RGEs which include the full effect of kinetic mixing 
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The running parameters are then used to calculate the tree level mass spectrum. However, 
it is well known that the one-loop corrections can be very important for particular particles 
and have to be taken into account. The best known example is the light MSSM Higgs 
boson which get shifted by up to 50% per-cent in case of heavy stops. Similar effects can 
be expected in the extended Higgs sector especially since these can be very light at tree- 
level. Similarly, the gauginos arising in an extended gauge sector can be potentially light 



and receive important corrections at one-loop 60] . To take these and all other possible 



effects into account we use a complete one-loop correction of the entire mass spectrum. Our 



procedure to calculate the one-loop masses is based on the method proposed in Ref. 68 



first, all running DR parameters are calculated at the SUSY scale and the SUSY masses 
at tree-level are derived. The EW vevs Vd and v u are afterwards re-calculated using the 
one-loop corrected Z mass and demanding 



(9bl9l + 9bl9r + 9 2 l9r) 



2 



Z Z ~ %| L + gj) [ d + J ( } 

in addition with the running value of tan /3. Note that 5m 2 z as well as all other self-energies 
include the corrections originated by all particles present in the mBLR. These calculations 
are performed in DR scheme and 't Hooft gauge. Also the complete dependence on the 
external momenta are taken into account. The re-calculated vevs are afterwards used to 
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solve the tree-level tadpole equations again and to re-calculate the tree-level mass spectrum 
as well as all vertices entering the one-loop corrections. Using these vertices and masses, the 
one-loop corrections 5ti to the tadpole equations are derived and we use as renormalization 
condition 

ti - SU = . (A13) 

These one-loop corrected tadpole equations are solved with respect to the same parameter 
as at tree level resulting in new parameters Bjp , [i^, B^tl respectively fj,^\ Bjp , m^R , 

2 (1) 

m£ R . The final step is to calculate all self-energies for different particles and to use those 
to get the one-loop corrected mass spectrum. 

1. Real scalars: for a real scalar <j), the one-loop corrections are included by calculating 



(38 



the real part of the poles of the corresponding propagator matrices 

Det[p?l-mJ ilL (p 2 )] =0, (A14) 

where 

ml lL (p 2 )=ml T -U4p 2 ). (A15) 

Equation (IA14I) has to be solved for each eigenvalue p 2 = m 2 which can be achieved 
in an iterative procedure. This has to be done also for charged scalars as well as the 
fermions. Note, m 2 - is the tree-level mass matrix but for the parameters fixed by the 
tadpole equations the one-loop corrected values are used. 

2. Complex scalars: for a complex scalar rj field we use at one-loop level 

m 2 1 'M) = m 2 T "-Tl 71 (p 2 ), (A16) 

While in case of sfermions rhjP agrees exactly with the tree-level mass matrix, for 
charged Higgs bosons /i^ and Bjp or and has to be used depending on the 
set of parameters the tadpole equations are solved for. 

3. Majorana fermions: the one-loop mass matrix of a Majorana x fermion is related 
to the tree-level mass matrix by 



mUpD = M t~\ 



+Mf (eS t (p?) + e£(p?)) 



(A17) 



where we have denoted the wave- function corrections by E^, TP L and the direct one- 
loop contribution to the mass by TP S . 

4. Dirac fermions: for a Dirac fermion \1/ one has to add the self-energies as 

M? L (p 2 ) = M* - S+(p 2 ) - ^%[p 2 )M* - M*E+(p 2 ). (A18) 



36 



Note, this procedure agrees with the method implemented in SPheno 3.1.10 to calculate 
the loop masses in the MSSM as well as with the code produced by SARAH 3 . . 39 or later. 
However, there are small differences to earlier versions of SPheno as well as other spectrum 
calculators: the MSSM equivalent of condition eq (1A12I) is often solved in an iterative way 
using the one-loop corrected parameters from the tadpole equations to calculate 5m z until 
m| + 5m 2 z has converged. In this context also /i^ and BjP are used in the vertices entering 
the one-loop corrections. However, these steps mix tree- and one-loop level and break 
therefore gauge invariance: when we tried this approach the relation between Goldstone and 
gauge bosons mass is violated. However, the numerical differences in case of the MSSM 
turned out to be rather small. 

As example we give the necessary formulae to calculate the one-loop corrections to the 
tadpole equations and the scalar Higgs masses in appendix IA 31 



3. 1-loop corrections of the Higgs sector 



As discussed in section IA 21 we have calculated the entire mass spectrum at one-loop. 
For that purpose it is necessary to calculate all possible 1-loop diagrams for the one- and 
two-point functions. As example we here give the corresponding expressions for the one- 
loop corrections of the tadpoles as well as the self-energy for the scalar Higgs fields. For all 
other self-energies we refer to the output of SARAH 10 . The results are expressed via Passarino 
Veltman integrals (68|. The basic integrals are 

d n q 1 



A (m) = lQn 2 Q 4 ~ n 
B (p,m 1 ,m 2 ) = 167r 2 Q 4 " n J 



i (27i) n q 2 — m 2 + is ' 
d n q 



i (27r) n 



is 



(q — p) 2 — m 2 + ie 



(A19) 
, (A20) 



with the renormalization scale Q. All the other, necessary functions can be expressed by Aq 
and Bq. For instance, 

1 



B 1 (p,m 1 ,m 2 ) 



2p 2 



Ao(m 2 ) - A (mi) + (p +m 1 - m 2 )B (p, m l7 m 2 ) 



(A21) 



and 



(A22) 
(A23) 



F (p, mi, m 2 ) =A (m 1 ) - 2A Q (m 2 ) - (2p 2 + 2m\ - ml)B (p, m u m 2 ) 
G Q {p, mi, m 2 ) =(p 2 -m\- m 2 2 )B (p, m 1 , m 2 ) - A (mi) - A (m 2 ) 

The numerical evalution of all loop-integrals is performed by SPheno. With this conventions 
we can write the one-loop tadpoles as 



-A (m z ^T a .^ z + -A (m 2 ZR ^T (7iyZ ^ ZR + 3A (m 2 w ^T ai 



w+w- 



10 To get the model files of the mBLR which is not yet part of the public version of SARAH please send a mail 
to the authors. 
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a=l 

2 3 

+ 4Vi () fm 2 _V i .+ -_m 2 _ + 12 V A (m 2 d )t l . m 2 

a=l a=l 

3 3 

+ 4 £ A (m 2 o ) r^mj, + 12 £ A (m 2 a ) r^m^ 

a=l a=l 
1 4 1 4 6 



a=l a=l a=l 

6 



a=l a=l 

7 9 

a=l a=l 
9 

+ 2^A (ml a y aitVatVa ml a (A24) 
0=1 

with <Tj = (urf, cr u , or, cr.R)f . denotes the vertex of the three particles x, y z, while 

F WX yz will De used for four-point interactions. For chiral couplings we use T L as coefficient 
of the left and T R as coefficient of the right polarization operator. For instance, T adiZ ,z is 
the coupling of a pure down-type Higgs to a Z boson while _„ - corresponds to the 
left-chiral part of the interaction of a i?-Higgs to a neutralino of the second generation. The 
expressions for all vertices can be obtained with SARAH. 
Using these conventions the self-energy for the scalar Higgs fields reads 



+ ~B Q (p 2 ,m z ,m z ^)Tl.^ z T au z R ,z + ~^ B *(f ,™\ R ,™z^Ki,z R ,zJ*i,z R ,z R 

+ 2 Bo { p2 ' m w-^ m w-)K j ,w+,w- r ^ i> w+,w- + 2A (m|)r CTi)CTiiX> z + 2A (m 2 Zi ^Y (Tt ^ z ^ Zh 

2 

+ 4A (m 2 w -y ai!(Ti!W+>w - -Y. A ^{ m2 H-) T o^H+,H- 

a=l 

2 2 
a=l 6=1 

- 2 V m~+ V B f V \ m 2 _ , m 2 _) m~- ( T L * _ + ._r R _ + ._ + T R * _ + „_r L . + 

Xa V \r ' Xa ' Xfc / Xb V °"i'Xa>X(, C^Xa^Xf, O^Xl.Xj, <?i,Xa ,X h J 

o=l 6=1 

+ y y G fp 2 ,m 2 _,m 2 _) (r L *. + ._r L _+ ._ +r i? *_ + ._r fl _ + 

/-^ /-^ U V Xa ' X 6 / V a j'Xa,X b °i,Xa,X b <?j,Xa,X b <Ti,Xa ,X b ) 



o=l 6=1 
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- 6 Y Y B ° m da > m l) m d b (^*,d a ,d b T aM + A 

a=l 6=1 

3 3 

+ 3 Y G ° ( p2 ' m l' m l) ( r a* ^A^iAA + ^AA^Aa) 



4b 



a=l 6=1 
3 3 



i m e a i m e b j m eb\^ aj,e a ,e b ^ ai,e a ,e b ^ aj ,e a ,e b ^ a-i,ea,e b 



a=l 6=1 

3 3 



+ WG n (Vm 2 m 2 Vr L * r L +r R * r fi 

~ Z-/ Z^ ^O^F 5 '" e „) '"et / ^ o-^e^ej/ (Ti,e a ,e b ~ 1 <jj,ea„e- b ai,e a ,e b 



a=l 6=1 

3 3 



6 m u a Yj A) (p ? m u a i m u b ^} m u b \^aj,ua,uj^ai,u a ,u b + ^ ^u,,^. « a m & 



a=l 6=1 

3 3 



a=l 6=1 

4 ^4 



o=l 
4 4 



a=l 6=1 
1 4 4 



2 2 2 

P , m h a ,™h b 



a=l 6=1 



a=l 



b 

^^o(ml a Jr CT . iCT . : g ;) g o - 3^A)(m?Jr CTi:CTi , GS , fia 

a=l a=l 



6 6 



a=l 6=1 



a=l 6=1 



r* ,- jT t* t + > > Bnip 2 , ml , m~ IT* -« s IV g* g t 

<jj,d* a ,d b a-i,d*,d b / , / , " \r i e a i e b J <jj,e*,e b <r M e a ,ej, 

+ 3 Y Y B ° (p 2 > m i' m fi*) r ^,fis,G6 r *i.fi!, 

a=l 6=1 
7 7 

- Vm,o Vso^.mlcm-olm^ofr^-o ~or R - ~ + r R *~ ~ r L -o V 

Z—/ X a u \r ' Xa' X^y *b \ aj,Xa>X b <ri,Xi,X b aj,Xa>X b °i,Xa,X b/ 

o=l 6=1 

1 7 7 

+ 2 G ° ( p2 ' m *"' ( r ?.x2,x? r lx°,x? + r ?,M r £,M, 

a=l 6=1 

9 9 9 

- A o (<)r ff „ ff „^,^ + Y Y B ° if> m l b ) T lj^ b T ^K,* b 

a=l 
9 9 

~ ^ m,/ -S Bo ( p2 ' m ^' 



a=l 6=1 



i ""v b J ""v b I aj,Va,y b ai,u a ,u b 1 o-j,v a ,v b cri,v a ,v b 



0=1 6=1 
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9 9 



I 1 V ( r? m 2 m 2 \(v L * V L -A- V R * V R 

a=l 6=1 

2 4 

+ 2 5Z r ^H/+,H- r ^^+,H- i71 o(p 2 ,m^,m 2 H/ _) +^r;. j7 ^ ob r CT . 7) ^ i ,F (p 2 ,mi oj) ,c 

6=1 6=1 
4 4 



2 2 2 

~ , tUa ,m 7 



(A25) 



4. RGEs 

The calculation of the reno rmalization group equations performed by SARAH is based on 
the generic expression of [1091 ] . In addition, the results of [59| are used to include the effect 
of kinetic mixing. 

The /3 functions for the parameters of a general superpotential written as 



w(4>) = l -^M 3 + -rrtfafafa 



(A26) 



can be e asily obta ined from the shown results for the anomalous dimensions by using the 
relations 
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(5f = Y*fo" 7p fc) 



J Y 



/i p( V } 



(A27) 
(A28) 



For the results of the other parameters as well as for the two-loop results which we skip here 
because of their length we suggest to use the function CalcRGEs [] of SARAH. 

a. Anomalous dimensions 



(i) 



7, 



(i) 

l 



7 W 

'H d 
<XR 



1 (l2 (Y}Y d + YjY u ) - (l8g 2 L + 32g 2 + g 2 BL + g 2 BLR ) 1 

-\ (2<?1 + 9bl + 9blr) 1 + YX + YX 
l - (2Tr(F e F e t) - 3g 2 + 6Tr (y d Y}) - g\ - g\ BL 
~ (2Tr(ny„t) - 3g\ + 6Tr (y u Y*) - g\ - g 2 RBL 



(A29) 
(A30) 
(A31) 
(A32) 



2g 2 R - 2g 2 RBL + 2VG9bl9rbl + 2V6g BLR g R - ?>g 2 BL - 3g 2 BLR + 4Tr (y s Y^ 



(A33) 
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7 £ = 4 ( " 2 { 9 R + ^) + 2 ^9bl9rbl + 2vW*<7fl - 3^| L - 3g 2 BLR J (A34) 
7? = Y}Y, (A35) 

7 * } = 12 ( 24F " ^ " ( 2v ^ sitoBi + ^blbSr + 32 ^ + + 6 9 2 rbl + <?Il + 9 2 blr) l) 

(A36) 

7^ = ^ (24F;Fj - ( - 2VQg BL g RBL - 2VQg BLR g R + 32^ + 6g 2 R + 6g 2 RBL + g% L + g 2 BLR ) l) 

(A37) 

7? } = t ( - (^(qr + g 2 RB ^j - 2\flg BL g RBL - 2\flg BLR g R + 3g 2 BL + 3#| M ) 1 

+ 4(20^ + 0?)) (A38) 

7 ^ = \ ( - ( 2 {& + ^) + 2 ^9bl9rbl + 2VEg BLR g R + 3^| L + 3^) 1 + 8F e *r e T ) 

(A39) 



6. Gauge Couplings 



^sl = 4 ( 27 ^bl - 2 ^>9bl9rbl + £bl (27g 2 BLR + 30^ Bi - v^blr^ 

+ fi'BL/? (30£f R - V(igBLR) tfftBZ,) (A40) 



7 (llgBLgBLRgRBL + 27g 2 BLR g R - 2y/6g BLR g 2 R + 30#| + 30g R g RBL 

VQgBLgRgRBL - V^BL^IbL ) ( A41 ) 



4 



A?bm = j(9 BL^ 9 BLR ~ >/GgR) + 9bl(^9R9RBL ~ V6g B LR.gRBL^ 

+ gBLR (27g 2 BLR - 2y/§g BLR g R + 30<?|) ) (A42) 
^oBi = 4 ( 27 9bl9rbl + 30^ BL + S^bl) + ^^l (27 g BLR g R 

- V6 (2g 2 RBL + 0^) ) ~ V6gBLR9R9RBLj (A43) 



PS = 9t (A44) 
/#> = -3^ s 3 (A45) 
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c. Gaugino Mass Parameters 



Pm bl = 2 ( 2 7g 2 BL M BL -g BL (- 27g BLR M BR + 2V6g RBL M BL + VEg R M E 

+ g R BL (sOg RBL M BL + 30g R M BR - \fcg BLR M BR ) ) ( A46) 

^ = - (27g 2 BLR M R + 30 to ( 5m M Bfl + g R M R ^j + p B£ (27<?blk - \fcg R ) M BR 

- V6g B LR (2g R M R + g RBL M BR ^ ) (A47) 

Z^mL = 4 (^7g 2 BL M BR + 27g 2 BLR M BR - V6g BLR (2g R M BR + g RBL (m bl + M R ) ) + 30 (^g 2 R M BR 
+ g 2 RBL M BR + (m bl + M R ) ) + (27# BLfl (m bl + Mr 



rblM br + g R ^M BL + M fi j j j j (A48) 

'm 2 - (A49) 
Bfl = -6g 2 M 3 (A50) 



[1] M. Malinsky, J. C. Romao and J. W. F. Valle, Phys. Rev. Lett. 95 (2005) 161801 
[arXiv:hep-ph/0506296] . 

[2] V. De Romeri, M. Hirsch, M. Malinsky, Phys. Rev. D84 053012 (2011); [arXiv: 1107.34121 
[hep-ph]]. 

[3] R. Mohapatra and J. Valle, Phys. Rev. D34, 1642 (1986). 

[4] E. K. Akhmedov, M. Lindner, E. Schnapka and J. W. F. Valle, Phys. Rev. D 53, 2752 (1996) 
[hep-ph/9509255| . 

[5] E. K. Akhmedov, M. Lindner, E. Schnapka and J. W. F. Valle, Phys. Lett. B 368, 270 (1996) 
[hep-ph/9507275| . 

[6] M. Hirsch, M. Malinsky, W. Porod, L. Reichert and F. Staub, JHEP 1202, 084 (2012) 
|arXiv:1110.3037l [hep-ph]]. 

[7] ATLAS Collaboration, G. Aad et al. , Phys. Lett. B710, 49 (2012). [aTXiv: 1202 .14081 

[8] CMS Collaboration, S. Chatrchyan et al. , (2012). irXiv: 1202 .14881 

[9] L. J. Hall, D. Pinner and J. T. Ruderman. larXiv:1112.2~703l [hep-ph], 
[10] H. Baer, V. Barger and A. Mustafayev, larXiv: 11 12.30171 [hep-ph]. 
[11] J. L. Feng, K. T. Matchev and D. Sanford. [arXiv:1112.302l1 [hep-ph], 
[12] S. Heinemeyer, O. Stal and G. Weiglein, la7Xiv:1112.3026l [hep-ph], 

[13] A. Arbey, M. Battaglia, A. Djouadi, F. Mahmoudi and J. Quevillon, Phys. Lett. B 708 
(2012) 162 [arXiv:1112.3028l [hep-ph]]. 



42 



[14] A. Arbey, M. Battaglia and F. Mahmoudi, Eur. Phys. J. C 72 (2012) 1906 |arXiv:1112.3032l 
[hep-ph]]. 

[15] P. Draper, P. Meade, M. Reece and D. Shih, arXiv:1112.3068l [hep-ph]. 

[16] T. Moroi, R. Sato and T. T. Yanagida, Phys. Lett. B 709 (2012) 218 [ arXiv:1112.3l42l 
[hep-ph]]. 

[17] M. Carena, S. Gori, N. R. Shah and C. E. M. Wagner. larXiv:1112.3336l [hep-ph]. 

[18] IT. Ellwanger. larXiv:1112.3548l [hep-ph], 

[19] O. Buchmueller, et al. larXiv:1112.3564l [hep-ph]. 

[20] S. Akula, B. Altunkaynak, D. Feldman, P. Nath and G. Peim. larXrv:1112.3645l [hep-ph]. 

[21] M. Kadastik, K. Kannike, A. Racioppi and M. Raidal. larXiv: 11 12 .36471 [hep-ph]. 

[22] J. Cao, Z. Heng, D. Li and J. M. Yang. larXw:1112.4391l [hep-ph], 

[23] A. Arvanitaki and G. Villadoro. larXiv: 11 12 .48351 [hep-ph]. 

[24] M. Gozdz. larXiv: 1201 .08751 [hep-ph], 

[25] J. F. Gunion, Y. Jiang and S. Kraml. larXiv:1201.0982l [hep-ph], 

[26] G. G. Ross and K. Schmidt-Hoberg. larXiv:1108.1284 [hep-ph]. 

[27] G. G. Ross, K. Schmidt-Hoberg and F. Staub. la7Xrv:1205.1509l [hep-ph]. 

[28] P. Fileviez Perez. larXiv:1201. 15011 [hep-ph], 

[29] N. Karagiannakis, G. Lazarides and C. Pallis. larXiv:1201.2111l [hep-ph]. 

[30] S. F. King, M. Muhlleitner and R. Nevzorov. larXiv:1201.26"7T1 [hep-ph]. 

[31] Z. Kang, J. Li and T. Li, larXiv:1201.5305l [hep-ph]. 

[32] C. -F. Chang, K. Cheung, Y. -C. Lin and T. -C. Yuan. larXiv: 1202.00541 [hep-ph], 

[33] L. Aparicio, D. G. Cerdeno and L. E. Ibanez. larXiv:1202.0822l [hep-ph]. 

[34] L. Roszkowski, E. M. Sessolo and Y. -L. S. Tsai, arXiv: 1202.15031 [hep-ph]. 

[35] J. Ellis and K. A. Olive. la7Xw:1202.3262l [hep-ph]. 

[36] H. Baer, V. Barger and A. Mustafavev. larXiv:1202.4038l [hep-ph]. 

[37] N. Desai, B. Mukhopadhyaya and S. Nivogi. larX"iv:1202.51901 [hep-ph]. 

[38] J. Cao, Z. Heng, J. M. Yang, Y. Zhang and J. Zhu. larXiv: 1202.5821 1 [hep-ph], 

[39] L. Maiani, A. D. Polosa and V. Riauer. la?Xiv:1202.5998l [hep-ph]. 

[40] T. Cheng, J. Li, T. Li, D. V. Nanopoulos and C. Tong. larXiv:1202.6088l [hep-ph]. 

[41] N. Christensen, T. Han and S. Su. larXiv:1203.3207l [hep-ph], 

[42] D. A. Vasquez, G. Belanger, C. Boehm, J. Da Silva, P. Richardson and C. Wymant, 

larXiv:1203.3446l [hep-ph]. 

[43] U. Ellwanger and C. Hugonie. larXiv: 1203.50481 [hep-ph], 

[44] Updated results can be found on the web-page: https:/ /twiki.cern.ch/twiki/bin/view/AtlasPublic; 

see, especially: ATLAS-CONF-2012-033 

[45] Updated results can be found on the web-page: 

https: / /twiki.cern.ch /twiki /bin / view / CMSPublic /PhysicsResults 

[46] H. E. Haber, M. Sher, Phys. Rev. D35 (1987) 2206. 

[47] M. Drees, Phys. Rev. D35 (1987) 2910-2913. 

[48] M. Cvetic, D. A. Demir, J. R. Espinosa, L. L. Everett, P. Langacker, Phys. Rev. D56 (1997) 



43 



2861. [hep-ph/9703317] . 

[49] Y. Zhang, H. An, X. d. Ji and R. N. Mohapatra, Phys. Rev. D 78, 011302 (2008) 
|arXiv:0804.0268l [hep-ph]]. 

E. Ma, [arXiv: 1108 .40291 [hep-ph]]. 

K. Nakamura et al. [Particle Data Group Collaboration], J. Phys. G G 37, 075021 (2010). 
T. Falk, K. A. Olive and M. Srednicki, Phys. Lett. B 339, 248 (1994) [hep-ph/9409270] . 
A. Ali, D. A. Demir, M. Frank and I. Turan, Phys. Rev. D 79 (2009) 095001 [arXiv:0902.3826l 
[hep-ph]]. 

A. Belyaev, J. P. Hall, S. F. King and P. Svantesson. larXiv:1203.2495l [hep-ph], 
S. P. Martin, Phys. Rev. D 46, 2769 (1992) [hep-ph/9207"2T8] . 

B. C. Allanach et al, Comput. Phys. Commun. 180 (2009) 8-25. [arXiv: 080 1.00451 [hep-ph]]. 
B. Holdom, Phys. Lett. B 166, 196 (1986). 

K. S. Babu, C. F. Kolda and J. March-Russell, Phys. Rev. D 57 (1998) 6788 
[arXiv:hep-ph/9710441~] . 

R. Fonseca, M. Malinsky, W. Porod, F. Staub, Nucl. Phys. B854 (2012) 28-53. 
|arXiv:1107.2670l [hep-ph]]. 

B. O'Leary, W. Porod and F. Staub, JHEP 1205 (2012) 042 [arXiv:1112.4600l [hep-ph]]. 

F. Braam and J. Reuter. larXiv:1107.2806l [hep-ph], 

E. J. Chun, J. -C. Park and S. Scopel, JHEP 1102 (2011) 100 [arXiv: 101 1.33001 [hep-ph]]. 
Y. Mambrini, JCAP 1107, 009 (2011) [arXiv: 1104.47991 [hep-ph]]. 
T. G. Rizzo, Phys. Rev. D 85, 055010 (2012) [arXiv: 120 1.28981 [hep-ph]]. 
D. Matalliotakis and H. P. Nilles, Nucl. Phys. B 435, 115 (1995) jhep-ph/9407251] . 
N. Polonsky and A. Pomarol, Phys. Rev. D 51, 6532 (1995) [hep-ph/941023l] . 
V. Berezinsky, A. Bottino, J. R. Ellis, N. Fornengo, G. Mignola and S. Scopel, Astropart. 
Phys. 5, 1 (1996) [hep-ph/gSOS^ . 

D. M. Pierce, J. A. Bagger, K. T. Matchev, R. -j. Zhang, Nucl. Phys. B491 (1997) 3-67. 
[hep-ph/9606211| . 

W. Porod, Comput. Phys. Commun. 153 (2003) 275 [arXiv:hep-ph/0301101] . 
W. Porod and F. Staub. larXiv:1104.1573l [hep-ph]. 

F. Staub. la7Xiv:0806.0538l [hep-ph], 
F. Staub. larXiv:0909.2863l [hep-ph], 

F. Staub, Comput. Phys. Commun. 182 (2011) 808 |arXiv:1002.0840 i [hep-ph]]. 
ATLAS conference note ATLAS-CONF-2012-007, see [44] 
J. Erler, P. Langacker, S. Munir, E. Rojas, [arXiv:1103.2659l [hep-ph]]. 

M. Hirsch, H. V. Klapdor-Kleingrothaus and S. G. Kovalenko, Phys. Lett. B 398, 311 (1997) 
[hep-ph/9701253| . 

M. Hirsch, J. W. F. Valle, W. Porod, J. C. Romao and A. Villanova del Moral, Phys. Rev. 
D 78 (2008) 013006 |arXiv:0804.4072l [hep-ph]]. 

M. Hirsch, S. Kaneko and W. Porod, Phys. Rev. D 78, 093004 (2008) [arXiv:0806.336T1 
[hep-ph]]. 



44 



[79] J. N. Esteves, J. C. Romao, M. Hirsch, F. Staub and W. Porod, Phys. Rev. D 83, 013003 



(2011) [arXiv: 1010.60001 [hep-ph]]. 

M. Hirsch, F. Staub and A. Vicente, larXiv: 1202.18251 [hep-ph]. 

J. Adam et al. [MEG collaboration], Phys. Rev. Lett. 107 (2011) 171801 [arXiv: 1107.55471 
[hep-ex]]. 

P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein and K. E. Williams, Comput. Phys. Com- 
mun. 181 (2010) 138 |arXiv:081 1.41691 [hep-ph]]. 

P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein and K. E. Williams, arXiv:1102.1898l [hep- 
ph]. 

G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 611 (2001) 403. 
A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 631 (2002) 195. 
A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 643 (2002) 79. 
A. Dedes and P. Slavich, Nucl. Phys. B 657 (2003) 333. 

A. Dedes, G. Degrassi and P. Slavich, Nucl. Phys. B 672 (2003) 144. 

B. C. Allanach, A. Djouadi, J. L. Kneur, W. Porod and P. Slavich, JHEP 0409 (2004) 044 
[hep-ph/04061661 - 

A. Nisati, for the ATLAS collaboration, talk presented at Lepton-Photon Conference 2011, 
Mumbai, India, August 2011. 

V. Sharma, for the CMS collaboration, talk presented at Lepton-Photon Conference 2011, 

Mumbai, India, August 2011. 

E. Etzion, private communication. 

G. Aad et al. [ATLAS Collaboration], Phys. Rev. Lett. 107 (2011) 272002 [arXiv: 1108. 15821 
[hep-ex]]. 

S. Chatrchyan et al. [CMS Collaboration], JHEP 1105 (2011) 093 |arXiv:1103.098T1 [hep-ex]]. 
W. Kilian, T. Ohl and J. Reuter, Eur. Phys. J. C 71, 1742 (2011) [ arXiv:0708.4233l [hep-ph]]. 



M. Moretti, T. Ohl and J. Reuter, arXiv:hep-ph/0102195 



F. Staub, T. Ohl, W. Porod, C. Speckner, [arXiv: 1109. 51471 [hep-ph]]. 

B. Fuks, M. Klasen, F. Ledroit, Q. Li and J. Morel, Nucl. Phys. B 797 (2008) 322 
|arXiv:071 1.07491 [hep-ph]]. 

F. del Aguila, J. A. Aguilar-Saavedra and R. Pittau, JHEP 0710 (2007) 047 
[hep-ph/0703261] . 

J. A. Aguilar-Saavedra, Nucl. Phys. B 828 (2010) 289 |arXiv:0905.222T1 [hep-ph]]. 
S. P. Das, F. F. Deppisch, O. Kittel and J. W. F. Valle. larXiv:1206.0256l [hep-ph], 
S. Gopalakrishna, A. de Gouvea and W. Porod, JCAP 0605 (2006) 005 [hep-ph/0602"027j . 

C. Arina and N. Fornengo, JHEP 0711 (2007) 029 [arXiv:0709.4477l [hep-ph]]. 

Z. Thomas, D. Tucker-Smith and N. Weiner, Phys. Rev. D 77 (2008) 115015 [arXiv:0712.4T46l 
[hep-ph]]. 

D. G. Cerdeno, C. Munoz and O. Seto, Phys. Rev. D 79 (2009) 023510 |arXiv:0807.3029l 
[hep-ph]]. 

P. Bandyopadhyay, E. J. Chun and J. -C. Park, JHEP 1106 (2011) 129 [arXiv: 1105. 16521 



45 



[hep-ph]]. 

[107] G. Belanger, J. Da Silva and A. Pukhov, JCAP 1112 (2011) 014 |arXiv:1110.24i~l [hep-ph]]. 
[108] B. Dumont, G. Belanger, S. Fichet, S. Kraml and T. Schwetz. la7Xiv:1206.152~T1 [hep-ph]. 
[109] S. P. Martin and M. T. Vaughn, Phys. Rev. D 50 (1994) 2282 [Erratum-ibid. D 78 (2008) 

039903] [arXiv:hep-ph/9311340| . 
[110] P. C. West, Phys. Lett. B 137 (1984) 371. 

[Ill] D. R. T. Jones and L. Mezincescu, Phys. Lett. B 138, 293 (1984). 



46 



